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"MICROSCOPIC STUDY OP THE ONSET OP LARGE DEFORMATIONS AND 
SHAPE- COEXISTENCE IN THE ZIRCONIUM REGION" 

It is generally felt that the neutron-proton (n-p) 
effective interactions possess deformation-producing tendency 
and the neutron- neutron (n-n) or proton-proton (p-p) effective 
interactions are mostly of the spherifying nature. These ideas 
have played pivotal roles in the development of the stretch 
scheme of Danos and Gillet / the roton model of Arirna and Gillet 
and the recent interacting boson model. 

The above assumptions regarding the nature of the two 
components of the effective interactions owe their origin 
mostly to the fact that the nuclei in the 2s— Id and 2p-lf shells: 
in which both protons and neutrons occupy half-filled sub- shells' 
usually display a number of features associated with well- i 

deformed intrinsic states. On the other hand, the low-lying j 
spectra in the neutron-excess nuclei in these shells are 

generally marked by an absence of deformation characteristics, j 

20 44 60 . ; 

Whereas the low-lying levels in Ne, Ti and Zn are quasi— j 

rotational, those in ^°0, 44 Ca and ^Ni display no analogous j 


deformation trends. 



Some time ago a new region of deformation has been 

discovered around A = 100; well-developed rotational spectra 

were observed in several highly neutron-rich isotopes of 

Zirconium and Molybdenum during a study of the fission fragments 
252 

of Cf. The observed B(E2) values were as enhanced as in the 
rare-earth and the actinide regions. 

It is difficult to ■understand the occurrence of rotational 
spectra in the Zirconium and Molybdenum isotopes in the frame- 
work of the traditional notions concerning the nature of the 
like-particle effective interactions. The recent observations 
in the Zirconium region thus call for a careful assessment of 
the deformation-producing tendency of the n-n (or p-p) channel 
of the effective interaction. 

A Hartree-Fock oriented analysis of the nature of the 
like-particle part of the two-body interaction is the starting 
point of the present work. The study provides a qualitative 
understanding of the occurrence of large deformations in the 
Zirconium region. This is followed by a microscopic study ofj 
number of interesting features associated with the observed 

| 

quasirotational spectra in the neutron-rich isotopes of f 

Zirconium, Molybdenum and Palladium. ' 

Chapter I presents the background and the genesis of th? 
topics that we have selected for investigation, in the present 
work. 

In Chapter 1 1 we have given an analysis of the n-p part 
and the part operative between the like-particles of the ; 



xiv 


various available effective interactions operative in 2s-ld, 
2p-lf and 2p-lf-lg shells, as regards their deformation produ- 
cing and pairing tendencies. Contrary to the widely-held 
belief, we found that like-particle interactions have almost 
the same deformation-producing character as the neutron-proton 
ones. Deformed Hartree-Fock method has been used for extracting 
the information concerning the deformation-producing tendency 
of the interactions. In this connection, we have examined the 
reliability of the various HF states involved by carrying out 
second-order perturbative calculations. A significant result 
that emerges from the present study is the following. It turns 
out that the relative irrportance of the pairing effects in the 
like-particle as well as the neutron-proton channels depends 
critically on the overall spread of the spherical single partic! 
energies. This interplay of the single-particle energies and 
the various components of the effective interactions is shown 
to explain the absence of rotational behaviour in the neutron- 
excess nuclei in the 2s- Id and 2p-lf shells. 

The above analysis of the various parts of the effective 
interactions proves instrumental in explaining the dramatic 
emergence of the rotational spectra in the A >100 region. 
Chapter III presents a systematic study of the onset of large 
deformations in the highly neutron-rich Zirconium and Molybdenur 
isotqpes. The calculations employ pairing-plus-quadrupole- j 

quadrupole effective interactions operating in a valence space | 

i 

1 
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XV 


80 

orbits for protons as well as neutrons. The nucleus Zr 

has been considered to be an inert core. The present choice 

of the large valence space ensures an unbiased study of the 

deformation effects in this region. The yrast spectra in 

various cases have been obtained in the framework of the 

Cranked Hartree-Fock-Bogoliubov approach. A number of 

interesting observed features in the Zirconium region connected 

with the dramatic onset of large deformations as well as the 

coexistence of deformed and spherical intrinsic shapes can be 

successfully understood in terms of the important roles played 

94 

by the deformation of Sr core and the increased participa- 
tion of the lh 1 ^2 level i n the valence space. 

A new phenomenon has some time ago been observed in the 
102 100 

level schemes of Pd and Pd populated by heavy-ion 
induced reactions r the ground state bands of these nuclei 
reach J = 14, but "fork" into two branches, at J = 8 + state 
in ^^Pd and at the J = 6 + state in "^^Pd. Further, the ground 
state bands in these nuclei have been observed to reveal 
significant departures from the rotational behaviour. 

Although attempts have been made recently to explain the 
observed spectra phenomenologically by invoking the "variable- 
moment of inertia" (VMI) model, microscopic calculations have 
been lacking. The successful calculation of the yrast spectra I 
in Zr and Mo isotopes encouraged us to make an attempt. at ; 

explaining the yrast and the yrare spectra in the nuclei 

100 102 ’ 
* Pd employing the 3s-'2d~lg-lh,, 1 /0 configuration space, , 



xvi 


in conjuction with the pairing-plus-quadrupole-guadrupole 
effective interaction. The results are discussed in the 
Chapter IV. Our first-ever microscopic study of Palladium 
isotopes brought out again the important role played by the 
lhn/2 orbital vis-a-vis the occurrence of forking effect as 
well as the significant deviations from quasirotational 
behaviour. Further, in sharp contrast the suggestions made 
in a couple of phenomenological studies of Pd isotopes in 
the recent past, it turned out that it is unnecessary to 
invoke asymmetric deformations in connection with the 
observed features of the yrast and the yrare spectra in the 
Palladium isotopes. 



CHAPTER I 


INTRODUCTION 

It is generally found that the nuclei in the 2s- Id and 

2p-lf shells, in which both protons and neutrons occupy half- 

filled sub-shells, display a number of features which are 

distinctly interpretable in terms of deformed intrinsic 

20 24 44 52 60 

states. Nuclei such as Ne, Mg, Ti, Fe and Zn are 
some typical examples. On the other hand, the nuclei consist- 
ing of only neutrons as valence particles display a conspicuous 

lack of rotational collectivity. Some examples of this category 

20 22 24 

are provided by the Oxygen isotopes ( * * 0), Calcium 

44—50 58—64 

isotopes ( Ca) and Nickel isotopesC Ni). 

An attempt is usually made to rationalize the above- 

mentioned features by associating the deformation-producing 

and the pairing tendencies with the neutron-proton ( ) and 

the like-particle ( v nn+ pp) effective interactions respectively, 

■In fact these assumptions concerning the two channels of 

effective interactions are implicit in the formulation of the 

1 

stretch scheme of Danos and Gillet as well as the roton 

2 

model of Arima and Gillet . The idea that the neutron-proton 

interaction is the major source of deformations is also inherent. 

3 

to some degree, in the recent Interacting Boson Model . 

4 

Recently Cheifetz et al. have discovered a new region of 
deformation around A=100. By studying the fission fragments 



2 


252 

of Cf, they observed rotational spectra in the nuclei 
l00/102 Zr x02 # 104, 106 Mq ^ The measured B(E2) values for 

these nuclei were as enhanced as in the rare-earth and actinide 
regions. 

90 

The assumption that the nucleus Zr can be considered 
"inert", has formed the basis of a large number of reasonably 
successful shell-model studies^ in the nuclei with A=91-96. 

The discovery of the we 11- developed rotational spectra in the 
highly neutron-rich Zr isotopes, therefore, appears quite 
puzzling? it seems very difficult to reconcile the observed 
onset of large deformations in these isotopes with the tradi- 
tional pairing caricature of the like-particle interactions. 

The recent experimental observations in the Zirconium 
region thus cast doubts on the correctness of the usual notions 
which associate the pairing or spherifying character with the 
like-particle effective interactions. It thus appears worth- 
while to carry out an assessment of the deformation-producing 
role of the like-particle part of the two-body effective 
interactions operative in the 2s- Id, 2p-lf and 2p-lf-lg shells. 
This is precisely what we have done in Chapter II, where the 
Hartree-Pock (HP) method has been employed to glean information 
about the deformation-producing roles of the interactions V 

and V . Our calculations differ from those carried out by 

nn+pp ■* 

6 

Nair et al. some time ago in that we have made an attempt at 
delineating the roles played by the spherical single-particle 



Hamiltonian. Our results bring out the link between the nature 
of the spectrum of spherical -SPE's and the possibility of an 
onset of rotational features. We find that the interaction 
V nn+pp ■ tlas a k° ut the same deformation-producing tendency as the 
interaction V and that they both have a dominant quadrupole- 
quadrupole character. We have also tried to examine, by 
carrying out second-order perturbative calculations, the 
connection between the overall spread of the SPE's and the 
relative importance of the pairing effects in the like-particle 

7 

as well as the neutron-proton channels. The calculations show 
that it is erroneous to interpret the absence of rotational 
features in the neutron- excess nuclei in the 2s-ld and 2p-lf 
shells in terms of the "pairing" nature of the interaction 

V nn+pp * 

The above analysis of the various parts of the effective 
interactions proves instrumental in understanding, at least 
qualitatively, the occurrence of large deformations in the 
nuclei ^^^Zr/ the significant degeneracy of the relevant 

spherical SPE's for the A=80-100 nuclei offers the interesting 
possibility of building up sizable deformations in this region 
for even those nuclei which possess only neutrons as the valence 
particles. 

Chapter III presents a quantitative study of the onset 
of large deformations in the highly neutron-rich Zirconium 
and Molybdenum isotopes. We have employed the Cranked Hartree- 
Pock-Bogoliubov (CHFB) theory in conjunction with the pairing- 



4 


lus-quadrupole-quadrupole effective interactions operating in 

l valence space spanned by the 2p^ 2 *■ 1 %/2 * ' ^ s l/2 ' 

'^3/2 ' ^l/"> an ^ 1 ^ 1 ll/2 or ki ts f° r protons as well as neutrons. 
80 

Dhe nucleus Zr has been considered to be an inert core. The 
present choice of the large valence space ensures an unbiased 
study of the deformation trends in this region. 

, 8-11 


Federman and Pittel 


have recently conjectured that the 


100 102 

dramatic onset of large deformations in the nuclei ' Zr 

may be related to the deformation-producing role of the 

[(lg^ 2 )p " ( lg^ 2 ^^ n ] interactions. Their suggestions imply 

a correlation between the occurrence of deformations and the 

simultaneous occupation of the orbits (Ig^y^n anc ^ ^ 1 %/2^p* 

12 

Our explicit HFB calculations , however, reveal that the 
interaction of the lg^ 2 -protons with the lg^ 2 -neutrons plays 


11/2 


very little role. We find that the involvement of the lh 

orbit is one of the major factors governing the intrinsic 

100 102 

quadrupole deformations in the nuclei 'Zr and 
102,104,106 Mo> 

Recently an attempt was made by Federman, Pittel and 
10 

Canpos to examine the co-existence of the spherical as well 

as deformed shapes in the Zirconium region by carrying out 

shell-model studies of some of the lighter Zr isotopes 

94 

(94<A<100). Their calculations involved an inert Sr core 
with the protons confined to the orbits (2p^ 2 , -^ 9/2 * 

2dj-/ 2 ) and the neutrons to the orbits (33^2 • 2d 3/2 ' lg 7/2 ' 
ih. Our CHFB calculations however reveal that the 



5 


ccurrence of deformed intrinsic shapes in the Zirconium and 

olybdenum isotopes with A>100 is not compatible with the 

94 

ssumption of an inert Sr core. The yrast spectra that we 
iave obtained for the nuclei 100 ' 10 ^z r an d o by 

employing the usual Cranked Hartree-Fock-Bogoliubov (CHFB) 

i 3 14 

:echnique ' ' are in reasonable agreement with the experiments. 

Scharff-Goldhaber et al^ 5 have some time ago reported on 

the level schemes of the nuclei 100 Pd (N=54, 2=46) and 102 Pd 

populated by heavy-ion induced reactions. An important feature 

that characterizes the yrast (as well as the yrare) spectra in 

these isotopes is the significant departure from rotational 

behaviour. Further, the ground state bands of these nuclei 

reach J=14, but "fork" into two branches, at J 11 =8 + in '^ 2 Pd 
100 

and at J=6 in Pd. 

Although a couple of phenomenological studies have been 

made some time ago to explain the observed spectra in 

100,102pch a lack of any explicit microscopic calculation has 

prevented a clear understanding of the observed features. The 

"variable moment of inertia" (VMI) model has been partially 

successful in explaining the observed deviations from the J(J+1) 

17 

law in the Palladium isotopes. Further, Smith and Volkov 

* 

have invoked the instability of the symmetric deformed nuclei 
towards asymmetric deformations in order to explain phenomeno- 
logically the forking effect. 

The successful qualitative study of the yrast spectra in 
the doubly even Zr and Mo isotopes encouraged us to make an 



6 


ttempt at providing a microscopic understanding of the yrast 

nd the yrare spectra in the Pd isotopes by carrying out the 

18 

ariation-after-pro jection (VAP) calculations with the HFB 
nsatz for the intrinsic wavef unction. Here again we employ 
.he pairing-plus-quadrupole-quadrupole effective interaction 
perating in the 3s-2d-lg-lh^^2 configuration space. Although 
>oth the CHFB and the VAP procedures allow for the possibility 
>f having different intrinsic state for each yrast level, the 
calculations of the intra-band electromagnetic transition 
>robabilities are easier to carry out if one employs the latter 
nethod. 


Our VAP results presented in Chapter IV for the nuclei 
100, 102p^ reveal a somewhat dichotomic character of the yrast 
spectra in these nuclei; we observe a sudden change of the 
intrinsic state in going from the yrast J =4 + to J =6 + . The 
calculation is successful in reproducing to a large extent, 
the observed deviations from the J( J+l) —law. 


The microscopic study of the Palladium isotopes that we 

present in Chapter IV brings out the irrportant role played 

by the lh^-j^ vis-a-vis the significant departure of the 

yrast band from the quasi rotational behaviour as well as the 

appearance of the yrare levels with J>6 at low excitation 

energies. Further, in contrast with the earlier phenomenolo- 
17 

gical studies , it turns out that it is not necessary to 


invoke asymmetric deformations in connection with the yrast 
and yrare spectra in Palladium isotopes. 



Finally, in • Chapter V we summarise the results 


btained in, and wisdom gained from Chapters II, III 

n a IV, 
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CHAPTER II 

DEFORMATION-PRODUCING TENDENCY OF THE 
LIKE-PARTICLE EFFECTIVE INTERACTIONS IN NUCLEI 

H*1 Introduction 

The nuclei in the 2s- Id and 2p-lf shells in which both 
protons and neutrons occupy half-filled sub-shells, usually 
display a number of features associated with we 11- deformed 
intrinsic states. On the other hand, the low-lying spectra in 
the neutron-excess nuclei in these shells are generally marlaed 
by an absence of deformation characteristics. Whereas the low- 
lying levels in 20 Ne, 44 Ti and 60 Zn are quasirotational, those 
in 2 ^0, 44 Ca and 6< ^Ni display no analogous deformation trends. 

As mentioned in the preceding chapter, the above observa- 
tions have led to a feeling that the neutron-proton (n-p) 
effective interactions possess deformation-producing tendency 
and the neutron-neutron (n-n) or proton-proton (p-p) effective 
interactions are mostly of the spherifying nature “ . These 

ideas have played a pivotal role in the development of the 

4 5 

stretch scheme of Danos and Gillet, the roton model of Arima 

and Gillet and the recent interacting Boson model^ of Arima 

et al . 

7 

Some time ago, Cheifetz et al. have observed well- developed 
rotational spectra in several highly neutron-rich Zr isotopes 
during a study of fission fragments of 


Cf. The observed 



10 


B(E2) values were as enhanced as in the rare-earth and actinide 
regions. It appears difficult to reconcile the observed dramatic 
onset of large deformations in the singly closed shell nuclei 
100 ' 1022 r w j_th the customary notion of "pairing" type n— n 
interactions. The observations in the Zr region thus call for 
a careful assessment of the deformation-producing tendency of 
the n-p as well as the n-n (or p-p) channels of the effective 
interactions. 

The deformation- or the field-producing character of a 
given two-body interaction (or a part of it) can quite unambi- 
guously be examined by carrying out a Hartree-Fock (HP) calcula- 
tion and by comparing the quadrupole moments of the resulting 
stable HF wave functions with the relevant S¥3 or the asymptotic 
Nilsson model predictions. This is precisely what we have done 
for some nuclei in the mass range A = 20-80 by employing 
various realistic effective interactions operative in the 2s-ld, 
2p-lf and 2p-lf-lg shells. 

An analysis of the deformation-producing tendencies of 

q g 

the T = 0 and T = 1 interactions by Sharma and Bhatt * had 
earlier revealed that these components have almost the same 
behaviour. The n-p interaction has, of course, both T = 0 and 
T = 1 parts. The like-particle interactions, on the other hand, 
constitute only a part of the full T = 1 interaction for 
systems with neutrons as well as protons. Therefore, it is not 
possible to glean information about the deforming nature of the 
n-n (or p-p) and n-p components from the earlier work. 
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In section II. 2 we give a brief outline of the deformed 

Hartree-Fock (HF) method. In section II. 3 we present the method 

for computing the two-particle two-hole (2p-2h) admixtures to 

the HF state in the framework of second-order perturbation 

theory. In section II. 4 we discuss the results obtained by 

carrying out deformed HF calculations by employing the different 

components of the two-body interaction. The significance of the 

conparison of the intrinsic quadrupole moments extracted from 

the two-body interactions (through the HF procedure) with the 

relevant SU3 values is pointed out. Our results show that^in 

contrast to the general belief, the like-particle conponent of 

the effective interaction possesses essentially the same 

def ormation-px'oducing character as the n-p component. Further, 

2 

a conparison of the computed<Q 0 > moments with the SU3 values shows 

2 2 

clearly the dominant quadrupole-quadrupole (q . q ) character 
of the various effective interactions that we have employed. 

Section II. 5 contains a discussion of the stability of 
the various HF minima involved in terms of the quantitative 
estimates of the 2p-2h admixtures. 

The results presented in section II. 6 bring out the 
important role played by the spectrum of the spherical singJLe- 
particle energies (SPE’s) vis-a-vis the relative importance of 
pairing effects in the like-particle and the neutron-proton 
channel. It turns out that the general absence of rotational 
features in the nuclei such as 2 ^0, ^Ca and ^Ni is not so 
much a reflection of the spherifying nature of the like-p article 
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interactions, it is just due to the effect of the significantly 
non-degenerate SPE ' s. 

It is easy to understand/ at least qualitatively/ the 
occurrence of large deformations in the Zirconium region in 
terms of the roles played by the spherical SPE's. This is 
discussed in section II. 7. Section II. 8 contains some conclud- 
ing remarks. 


II. 2 The Deformed Hartree-Fock Method . 

Theoretical discussions of the deformations in nuclei 

10 

are usually based on the unified model of Bohr and Mottelson 
in which single particle and collective rotational motions are 
distinguished and uncoupled. In this model/ e.jg. , the state 
of spin J belonging to the ground-state rotational band of an 
even-even nucleus is described by the product wave function 


D MO (a) IV 


The rotational motion is described by the symmetric top wave 
function 



(n) 


2J+1 


M 


(»Y) 


( 11 , 2 . 2 ) 


The single-particle motion is- described by the intrinsic state j 
|$ >. The wave function of the nucleus is only expected to ! 

factorize when the particle excitation energy or the spacing 

of intrinsic states is larger than the energy of collective j 

I 

excitation or level spacing in a rotational band. The conditions! 


I 
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are sometimes met in even-even N=Z nuclei in the A=20-60 mass 
region. 

The sirtplest assumption that can be made about the intrinsic 
state | $ > is that it is a Slater determinant of single nucleon 
orbits in a deformed field, 

|d>> = bt bt ... bt |o> (II. 2. 3) 

A 1 A 2 a A 

The notation of second quantization is used here/ and b^ is the 
fermion operator which creates a particle in the orbit X . It 
obeys the anti- commutation rules 

£ b X' b u3+ = 6 XyJ?V b yl+ = 0 (II. 2. 4) 

The purpose of the HF theory is to determine the wave 
functions of the orbits \ of the Slater determinant > . (II. 2. 3). 
This is achieved by requiring the energy of the system described 
by the wave function in eq. (II. 2. 3) to be stationary and 
minimum for infinitesimal variations of the orbits x . The 
energy of the system is obtained by taking, the expectation value 
of the Hamiltonian H in the state of eq* (II. 2. 3). The 
Hamiltonian H is defined 

A 2 

H = I V 1 Z V(i j) (II. 2. 5) 

1=1 2m 2 i/£j 

E = < <!> |H | $ > 

A A 

= I <x|t|x>+i Z < X y |V | X u > (II. 2. 6) 
X=1 Xy =1 

where is called the HF energy; 
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<x|t|x> i s the kinetic energy of the orbit X, 
and 

<Xu| v |x V >=~<Xy| V I y X > (II. 2 . 7 ) 

is an anti symmetri zed matrix element of the two-body interaction 
V(i j ) . 

It is advantageous to expand the orbits on a basis | jm > of 
known wave functions 


X > 


E 

j 


C x I jm > 
jm 1 


( 11 . 2 . 8 ) 


The basis j may be e.g. , the set |nljmT>of shell model states 
in a harmonic oscillator. Any basis may be. used to expand the 
orbits, provided one is able to calculate matrix elements of the 
interaction v with the basis states. The set of orbits X (both 
filled and empty) is assumed to form an orthonormal set of 
wave functions 


£ 

j 





E 

X 


* 


c y 

jm 


c 


X 

j 'm 1 



(11.2.9) 


The wave function jct>> is determined by the coefficients 
which become variational parameters. The energy will be 
stationary with the normalization condition (II. 2. 9) when 


^[<»IhU>- e x i c £ cjV]= o (II. 2. 10) 

jm 

Here e^ is introduced as a Lagrange multiplier. With the help 
of expansion (II. 2.8) of the orbits, the energy < <f> | H J > in 
eq. (II. 2. 10) may be expressed in terms of the known matrix 
elements of t and V in the known basis jjm>. The derivative 



(II. 2. 10) may then be calculated directly, and one obtains 

r A _». 

S L <jm|t|j , m'> + E . £ I C^' 

j'm» 


X=l j 2 m 7 J 1 m 1 <jm, j 1 m 1 |v( j'm' , J 2 m 2 > 


x C* ] c\ , = e, C . X 
J 2 m 2 J 'm* X jm 

oqs. (II. 2. 11) have the form of an eigenvalue problem 

I < jmlh Ij'mS c\ . = e, C * 

j« m . 1 |J J'm* X jm 

where h is a Hamiltonian given by its matrix elements 

A 

<jm|h|j , m'> = <jm[t|j , m’> + z < jm,x| V |j ‘m 1 , \ > 

X=1 


(II. 2. 11) 


(II. 2. 12) 


A 


X=1 j 1 m 1 j 2 m 2 J l m l 


= < jm |t |j 'm* > + E i E 

X=1 j ± r 

<jm, j^mJ vj j *m* , j 2 m 2 > C 


C * 


* 


X 

^ 2 m 2 


(II. 2. 13) 


where h is the HF single particle Hamiltonian. Equations 

(II. 2. 12) and (II. 2. 13) are the HF equations for the orbits X. 

11 

They may be solved by the following iteration procedures 

(a) An initial set of coefficients C ^ is guessed, 

jm 

and also which orbits are occupied. 


(b) With this set of coefficients the matrix elements 

of eq. (II. 2. 13) of the HF Hamiltonian are calculated. 

(c) Equation (II. 2. 12) is solved by diagonalizing the 
Hamiltonian matrix (II. 2. 13). A new set of coeffi- 
cients C.* is obtained, and A orbits are chosen to be 

jm 

filled by nucleons. One returns to step (b). This 
process is repeated until successive diagonalizations 
produce the same set of coefficients . 
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The iteration procedure is some^what complicated by the 
fact that there are many local minima of the energy (II. 2.6) and 
therefore many solutions of the HP eraiations. Various minima 
may be reached by various initial guesses, which therefore 
become crucial. Furthermore, after each diagonalization of the 
Hamiltonian matrix (II. 2. 13) one must decide which orbits to fill 
and it will be seen that the energies of the orbits. vary from 
one solution to another. 

In our calculations the starting wave functions have been 
generated froia the deformed single-particle Hamiltonian 

h = ej S.., + (II. 2. 14) 

We have tried to reduce the time involved in the calculations by 
storing all possible uncoupled matrix elements of the type 
< - ) "l m l / ^2 rn 2 I V { J 3 m 3 v J 4 m 4 > w 'bich are compatible with the axial 
symmetry of the orbits |x> . 

1 1 . 3 . Ca lculation of the Second-Order Admixtures to the HP state. 

The lowest-order correction to the HP intrinsic state 
will be through the admixture of 2p-2h type excited 
configurations to it. The 2p-2h configurations are limited to 
the HF single particle states of a given major shell. For the 
intrinsic HF state ] cf> > a typical 1 2p~2h > configuration is 
given by 

I 2p-2h > = b+ x b^ b x, b X2 |* > (II. 3. 1) 

The matrix element connecting the Hamiltonian (II. 2.5) 
with a 2p-2h configuration is then given by 



17 


<2p~2h|H|$ > = <y 1 y 2 |v|X 1 X 2 > (II. 3. 2) 

Here and y^ denote HP errpty orbits belonging to the same 
major shell as the occupied orbits X^ and X 2# 

The second order perturbation theory yields the following 
corrections for the HF-energy and the mass quadrupole moment. 

A = _ S S A „ < y- 1 ^ O l V l > / XT -3 ,\ 

y^<y 2 X-^<X 2 ^ 12 (Ix.3. 3) 


A Q (1 + C 2 p- 2 h ) 1 y 1 <V^ (A ^lV 2 x l x 2 )2 


x (q 


VI VI + q y 2 y 2 “ q X.X, “ q X 0 X. 


) 


l ,l l A 2 2 

The total wave function can be written as 


^ " C 2 p- 2 h^ ^HF + C 2 p- 2 h ^ 2 p- 2 h 


where we have 


C 1 

= 2 p- 2 h 

2 p- 2 h /, . - 1 2 

^ U+C 2 p- 2 h) 

pi _/S E a ) d 

2 P“ 2 h " V y 2 X X <X 2 -y lV2 X x ^ 


The amplitude A , . is given by 

y i V 2 12 

A X X = <U l U 2i V ' X l X 2 > A 

y l y 2 1 2 Ce !” +e -e . -e. ) 

y l y 2 X 1 X 2 


(II. 3.4) 


(II. 3.5) 


(II. 3.6) 


(II. 3.7) 


Here ^ y 2 1 V ^ X 1 X 2 > A an an ‘ t: l s y rru ’ ne1::r ^ zec l matrix element. In 
eq. (II. 3.4) the cross terms of the quadrupole operator between 
various 2p-2h configurations built on the HF-state have been 
neglected as their effect is expected to be very small. The 
matrix element of the quadrupole operator between the single 
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particle states X and is denoted by . 

1 2 

The squared amplitude = <*l* 2 p. 2h > 2 expresses 

the total intensity of the 2p-2h admixtures to the HP intrinsic 
states. 


I I • 4 Deformation-producing Tendency of the Like-particle 

Effective Interactions In The 2s- Id, 2p-lf and 
2p-lf-lg shells 

In order to assess the deformation or the field 

producing character of the parts V__ and V „ of various 
^ ^ nn nn+pp 

effective interactions operative in the 2s-ld, 2p-lf and 
2p-lf-lg shells, we have carried out deformed HF calculations 
and compared the intrinsic quadrupole moments of the stable HF 
minima with the relevant SU3 or the asymptotic Nilsson model 
predictions. 


In Table II. 1(a) we first present the results obtained in 

20 

the case of the nucleus Ne. The valence space is spanned by 

16 

the 1<%^2' 2s l/2 an< ^ 1< ^3/2 or kits. The nucleus 0 is treated 

as an inert core. The relevant effective two-body interaction 

12 

that we have employed is a renormalized G matrix due to Kuo 


which is a sum of G k are an< 3- G 3p_ih 


16 


0 core. The Kuo 
13 


interaction has earlier been errployed by Halbert et al. in 
a large number of successful shell-model calculations in the 
2s-ld shell. The one-body part of the Hamiltonian is specified 
by the observed 0 energies (in MeV) : e( 2s^) = 0.87, 
e(ld 5 y 2 ) =0.0 and e (13y 2 ) = 5.08. We coirpare here the 
quadrupole moments of the HF states resulting from 
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TABLE II. 1 

2 

The HP energies (E ) , intrinsic quadrupole moments <Q > and the 

20 ^ 

energy gaps (G) of the HP states of Ne obtained with different 

components (V) of the effective interaction. The quadrupole 

2 

moments are given in units of b , where b is the harmonic 

oscillator parameter. The SUG value for the cuadrupole moment 
20 2 3 

of Ne is 16.0 b . The last two columns summarize the results 
of the second-order perturbation calculation. 


V 

E(MeV) 

A 

o 

O to 

V 

G(MeV) 

&E 2p-2h 




(a) Experimental SPE 

1 s 


Pull 

-37.80 

15.45 

7.25 

-0.18 

00.95 

V 

nn+pp 

-20.45 

13. 76 

1. 15 

-1.08 

20.37 

V np 

-31.20 

15. 16 

4.89 

-0.66 

04.76 



(b) Degenerate SPE' 

s 


Pull 

-26.57 

15.58 

8.44 

-0. 15 

00.79 

V 

nn-i-pp 

-7.26 

15.61 

2. 31 

-0.49 

08.44 

£ 

-19. 30 

15.60 

6. 12 

-0.74 

05.36 



(i) the full Kuo interaction 

(ii) only the like-particle (n-n and p-p) part (labelled 
as v nn+ pp) °f "the interaction, with the n-p part 
set equal to zero, and 

(iii) only the n-p part (labelled as V ) of the interaction 

np 

with the like-particle part set equal to zero. 

It may be noted here that the calculations (ii) and (iii) 
violate isospin invariance. The main aim of our model calcula- 
tions is to isolate the deformation-producing tendency of the 

V and V parts of the effective interactions. The 

nn+pp np * 

results are not intended for comparison with any real data. 

Hence the violation of a well-established symmetry is not 
objectionable. 

2 

It is seen that the quadrupole deformation < Q Q > of the HF 

2 

state resulting from the full interaction is 15.45 b , where 

b ( = »Vm“) is the harmonic oscillator parameter. When only 

the like-particle interactions are employed, the quadrupole 

2 

d •-■formation is 13.76 b , whereas n-p interaction alone yields 
2 2 

< Qq >= 15.16 b . We thus find here that the n-p interactions 
yield only slightly more quadrupole deformation than the like- 
particle ones. 

With a given spectrum of the spherical SPE‘s the quadrupole 
deformations of the HF wave function depends not only on the 
multipole character of the interaction but also on its strength . 
We therefore switch off the SPE 1 s because our primary aim here 
is to isolate and compare the deformation-producing tendencies 
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of the n-p and the like-particle interactions. With the one- 
body part of the Hamiltonian absent, the wave functions (and 
hence the deformations) become indep endent of the strength of 
the interaction. The HP energy then depends linearly on the 
interaction strength. 

The results obtained with the degenerate SPE's (Table 

11.1(b)) reveal quite dramatically the identical deforming 

tendencies of the n-p and the like-particle interactions; 

2 2 

whereas the full interaction yields < Q Q > = 15.58 b , the 

intrinsic quadrupole- moments of the HP states resulting from 

the like-particle interactions and the n-p interactions are 
2 2 

15,61 b and 15,60 b respectively. 

Note also that the HP energy resulting from the full 
interaction is almost the sum of the energies obtained i n the 
case of the interactions v nn+ pp and V^. This again implies 
that the deformed HP wave functions obtained from the three 
interactions (viz. V nn # and the full interaction) are 

nearly identical. 

2 

The magnitude of the intrinsic <Q q > moments resulting from 

V and V are quite close to the SU3 values. To appreciate 

nn+pp np 1 

the significance of this result we first note that the HP 
self-consistent single-particle field obtained from a two-body 
interaction V can be written as 

V =Z 0 (E V ~ < a^ a * >}a* a« . 

H-P aBySQYpo Y o J a p 

(Here the expectation values <...•> are taken with respect to 
the HP state) , 



Assuming the potential V = V £ q(i).q(j) and neglecting 

i<j 

the exchange term (the Hartree approximation) we get 


V H = { V o J 2 <a| ( 3 u lB> < Ylq u |5> <a Y a 6 >}a a. a S* 

p y o ^ 


a 0 


Since q = E <a |q I e>at a 0 , V u can be recast in the form 
y a 0 1 a p H 

V„ = V Z (-l) M <q > q 

H o y ^y 

If we assume 14 that an axially symmetric, SU3 intrinsic 
state <>(Xo) of maximum weight has been employed in the r.h. s. 
of the above equation, then becomes proportional to q Q . 

Thus, within the Hartree approximation, the self-consistent 
intrinsic states resulting from a pure q. q two— body interaction 
is just the SU3 state of maximum weight. In view of this 
result, the closeness of our corrputed < Q q > values to the relevant 
SU3 values clearly implies the dominant q. q character of the 
three sets of interactions erployed, viz., the complete Kuo 
interaction, its n-p part and the part operative between the 
like particles. 


The validity of the above conclusions depends very much 
on the reliability of the HP wave functions involved. We shall 
return to the problem of assessing the goodness of the HP states 
later in the section II. 5. 


Next we present in Table II. 2 and II. 3 the results obtained 
for some nuclei in the 2p— If and 2p— If— Ig shells. The calcula- 
tions for the nuclei 44/48 TL, 48,52 Cr and 56 ' 62 Ni have been 

carried out by using the renormalized interaction, (G + G.^.^), 

15 

obtained by Kuo and Brown for the valence space consisting of 
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TABLE II. 2 

The HF energies as wall as the intrinsic quadrupole moments 
for some nuclei in the 2p-lf shell. The spherical SPE's have 
been taken to be degenerate. 







n 


Nucleus 

Interaction 

E(in MeV) 

< Q J > 

<Q ^SU3 


44 m . 

22 Tl 22 

Pull 

-16.42 

23. 24 


(a) 

V nn +V pp 

-4.64 

23. 28 

24.00 



V np 

-11.77 

23. 30 



48 Ti 

22 Xi 26 

Pull 

-30. 99 

35.37 


(b) 

V nn +V pp 

-11.03 

35. 39 

36.00 



V np 

-20.00 

35. 23 




Full 

-35.84 

35.32 


(c) 

48 Cr 

V +v 

-9.93 

35.09 

36.00 

24 24 

nn pp 






V np 

-25.95 

35.64 



52 Cr 

24 or 28 

Full 

-57.34 

41.57 


(d) 

V nn +V pp 

-16.33 

41.47 

42.00 



V np 

-41.09 

41.46 



56„. 

Full 

-94.68 

47.79 


(e) 

V +V 

-22.73 

47.31 

48.00 

2L 28 

nn pp 






V np 

-71.95 

47.79 



62 kt : 

28 Nl 34 

Full 

-127.08 

41. 92 


(f) 

V +V 
nn pp 

-29.28 

41.82 

42.00 



V 

-98.41 

41.02 




np 






TABLE II. 3- 
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The HF energies as well as the intrinsic quadrupole moments 
for some nuclei in the 2p-lf-lg shell. All the spherical 
SPE's are degenerate. 



Nucleus 

Interaction 

E 

<Q o > 

<Q o> max 



Pull 

-9.34 

20.13 


(a) 

60„ 

30 Zn 30 

V nn +V pp 

-2.42 

20. 14 

20.24 



V np 

-6.92 

20.14 




Full 

-26.87 

32.01 


(b) 

66 

30 Zn 36 

V nn +V pp 

-10.13 

29.91 

32.36 



V np 

-17.89 

31.84 




Full 

-26.99 

29.69 


(c) 

64 

3 2 Ge 3 2 

V nn +V pp 

-7.70 

29.23 

30.88 



V np 

-19.33 

29.94 




Full 

-61.26 

34.04 


(d) 

74 

32^42 

V nn +V pp 

-19. 90 

37.43 

38.48 



V np 

-44.17 

33.74 




Full 

-47.34 

36.98 


(e) 

68 Se 

34 &e 34 

V +V 
nn pp 

-13. 10 

36.32 

38.88 



V 

np 

-34.30 

37.34 




Full 

-82.54 

39.12 


(f) 

76 Se 

34 toe 42 

V nn +V pp 

-22.61 

40.96 

42.48. 



Vn P 

-61.65 

38.90 
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40 

the orbits 2p ^ 2 ,2P 1/2 an< ^ ^^5/2* T ^ e nuc l eus Ca I s 

considered as an inert core and the one-body part of the 

41 

Hamiltonian is defined by the observed Ca spectrum: 

(in MeV) e(if 7/2 )=o.O, e ( 2p 3/ , 2 ) =2. 1, e ( 2p 1 ^ 2 )=3„9 and 

60 66 

s(lf^ 2 ) = 6.5. The calculations for the nuclei * Zn, 
64/74 Ge and 68/76 Se have been carried out by choosing the 


valence space spanned by the 2p 


3/2' 4 ^5/2 / ^ p l/2 an< ^ 


9/2 


orbits. The interaction that we use is a sum of G, , G 0 

bare 3p-lh 

16 5 6 t 

and G 2p-2h ca l cu l ate( 3- fry Kuo x for the Ni core. The experi- 

57 

mental SPE’s, taken from Ni spectrum, are (in MeV): 
e(2py 2 ) = 0*0/ e ^ 2 Pi/2^ = l* 08 ' e ^^5/2^ ~ °*^8 anc ^ 
ff(l0gy^2^ ~ 2 »0* 


The minimum-energy HP states for the nuclei considered 

here possess axial symmetry. The striking similarity of the 

deformation-producing tendencies of the like-particle and the 

n-p interactions again manifests itself here in the (almost) 

2 

identical < Q Q > values resulting from these interactions for ■ 
various nuclei. Further, the HF quadrupole moments obtained 
from the n-p as well as the like-particle interactions are 
again almost the maximum possible ones for all the nuclei 
considered here, reflecting the dominant q, q character of the 
two components of the effective interaction. 


II. 5 Stability Of The Hartree-Fock Solutions 

The stability of the HF minima depends in a critical 
manner on the corrpetition between the magnitude of the energy 
gap (G) at the Fermi surface and the strength of the interaction. 



Usually an energy gap of about 2 MeV or more can be regarded 
safe enough to check the onset of pairing correlations in the 
HFB f rameAvork, in the 2s-ld or 2p— If shell. From the results 
presented in column 4, Table II. 1, one notes that this is indeed 
so, except in one case when we employ ( v nn , ) it the presence 

'ir ir 

of the experimental single-particle energies. 

We have also examined quantitatively the goodness of the 

HF states by computing all possible 2p-2h admixtures to them, 

in the framework of second-order perturbation theory. It should 

be mentioned here that the Kuo interaction employed in the 

present work has already been partially renormalized by 0 

core polarization processes. These processes cannot be 

described within the 2s-ld space. The 2p-2h states considered 

here correspond to the excitations from the occupied 2s-ld shell 

HF orbitals to the unoccupied 2s-ld shell orbitals. In columns 

5 and 6, Table II. 1, we give the lowering of the HF energies 

(AE), as well as the total intensity, C 0 ,,, (= I <tp j >2., 

zp- . 2p2h 2p2h J ' 

of the various 2p-2h admixtures. 

We first consider here the results obtained in the case (b) 

when the SPE’s are degenerate. In view of the smallness of the 

corrections obtained here, the HF wave functions obtained in 

20 

various cases for the nucleus Ne can be expected to provide 
a reasonably good approximation to the relevant exact shell * 

model wave functions. 

From the results obtained in the presence of the experimental 
SPE's (case a), we find that whereas the admixtures resulting 
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from v nn+ pp are about 20.4 percent, those due to V are only 
about 4.8 percent. As we shall see in the next section, this 
result is not confined to the present case alone; it is 
connected with the role of the spherical SPE's in giving rise 
to enhanced pairing effects in the neutron-excess nuclei in 
the 2s- Id, 2p~lf and 2p-lf-lg shells. 

From the results of the perturbation calculation for 
some of the axially symmetric nuclei in the 2p-lf and 2p-lf-lg 
shells (see columns 9-12, Table 11.4(a)) we note that various 
HF states obtained here can also be considered fairly reliable; 
the total intensity of all possible 2p-2h corrections is 
always less than 21 percent. It may be mentioned that the 
usual 2p-2h correlations involved in an HFB type calculation 
are only a small fraction of the total admixtures considered 
here - whereas one excites only time-reversed, K=0 pairs in 
the HFB scheme, we have included both K=0 as well as K/0 pairs. 

II . 6 Role of the Spherical SPE*s vis-a-vis the Enhanced 

Pairing Effects in the Like-particle Channel ™ 

We now turn to the question of the general absence 

of rotational collectivity in the neutron-excess nu’clei in 

A=l6-60 ■ region. In Table 11.4(b) we present the results 

44 48 60 

obtained in the case of the nuclei ' Ti and Zn by 
employing the relevant SPE's for the 2p-lf and 2p-lf-lg shells. 
A comparison of these results with those presented in 
Table 11.4(a) brings out dramatically the role of the experi- 
mental SPE's in these shells vis-a-vis the onset of pairing 
correlations in the neutron-excess nuclei. Whereas the 2p-2h 



I?he HF as well as the perturbation tneory uuudxucu m 

Interactions for some nuclei in the 2p-lf and 2p-lf-lg shells 
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admixtures for Ti resulting from the interaction V are 

nn+pp 

only 14 # 5 percent with the degenerate SPE's, they increase to 

57.5 percent v-hen the experimental SPE's are used. The 2p-2h 

admixtures resulting from V for the nuclei ^®Ti and ^Zn 

nn+pp 

are also seen to increase significantly in going from a degener- 
ate set of SPE 1 s to the experimental one. In sharp contrast to 
these results/ the 2p-2h admixtures due to V„ in fact decrease 
when we employ the experimental SPE 1 s in place of the degenerate 
ones. 

A number of Hartree-Fock-3ogoliubov calculations have been 
carried out in the recent past for the 2s- Id and 2p-lf shell 
nuclei. An important feature of these calculations is the 
treatment of pairing as well as the deformation effects on the 
same footing. These calculations have clearly established a 
link between the occurrence of non-rotational features and the 
extent of pairing correlations/ as reflected in the HFB occupa- 
tion probabilities and pairing gaps etc , in the N>Z nuclei in 
the 2s-ld and 2p-lf shells. 

The results presented here/ therefore, indicate that it 
is reasonable to attribute the observed departures from rotational 
behaviour in the neutron- excess nuclei in the A=l6-60 mass region, 
to the significantly non- degenerate spectrum of the SPE's in 
this region. 

II. 7 Occurrence of Large Deformations in the Zirconium Region 
The spherical SPE's in the Zr region are characterized 
by near- degeneracy; the spectrum of SPE's errployed in a number 
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of successful shell-model studies in this region by Vergados and 

17 18 

Kuo as well as by Federman and Pittel shows that five sub- 
shells, with a total occupancy of 32 like-particles, are clustered 
within just 4 MeV. In what follows we show that the emergence 
of the quasirotational spectra in this region can be understood, 
at least qualitatively, in terms of this feature of the SPE's. 

-4 

We have verified that the interaction v n n/pp = -0.02 b q. q, 

41 

when employed in conjunction with the Ca SPE's, can simulate 

2 

the KB interaction reasonably well as far as the <Q q > values for 

various Ca isotopes are concerned. The intrinsic HF states 

resulting from this interaction (as well as those resulting from 

the KB interaction) are near- spherical. We have employed this 

interaction in a valence space spanned by the lgg / / 2 / 2d 5/2 / 

2s 1/2 / 2< ^3/2' ^ 7/2 an< ^ ^11/2 or ^ ts f° r neutrons. The nucleus 
80 

Zr has been considered to be an inert core. The SPE's 

17 18 

employed are exactly the ones used in the earlier works ' , 

except for a slight lowering of the ener gY b Y MeV « 

The results of our HF calculations yield significantly small 

2 90—98 

< Qq > values for the isotopes "* Zr; the quadrupole moments 

2 

are always less than 18.7 b whereas the maximum possible value 

2 2 

for these isotopes ranges from about 47.6 b to 62.7 b for the 

given space. This is consistent with the observed near- spherical. 

100 102 

nature of these isotopes. For the nuclei * Zr, however, 

2 

we obtain significantly enhanced deformations; the computed < q q > 

2 

values are 53.3 and 64.7 percent of <Q n > rn _ >v respectively. 

vj max 
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II. 8 Conclusions 

To summarize, we have shown that the like-particle 
interactions, often characterized as the "pairing' 1 interactions, 
possess essentially the same deformation-producing tendency as 
the n-p interaction and that they both have a dominantly q. q 
character. Our calculations have revealed that the spectrum of 
spherical SPE's plays a very significant role vis-af-vis the 
relative importance of pairing correlations in the like-particle 
and the n-p channels. An appreciation of this role of the 
spherical SPE 1 s is found instrumental in reconciling the near- 
spherical nature of the N>Z nuclei in the 2s- Id and 2p-lf shells 
with the onset of large deformations in highly neutron-rich Zr 
isotopes. 
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* CHAPTER III 

MICROSCOPIC STUDY OP THE ONSET OF LARGE 
DEFORMATIONS IN THE ZIRCONIUM REGION 

III.1 Introduction 

In the last chapter we have mentioned that the 
significant degeneracy of the relevant spherical SPE's for 
A=80-100 nuclei offers the interesting possibility of build- 
ing up sizable deformations in this region for even those 
nuclei which have only neutrons as valence particles. In 
this chapter we present explicit calculations of the yrast 
spectra in the isotopes of Zirconium and Molybdenum. This 
has been done with a view to achieve, at a quantitative 
level, a microscopic understanding of the fascinating shape- 
transition that appears to be occurring in this region. 

In figure 1 we have given the observed low-lying spectra 
in some doubly-even isotopes of Zirconium and Molybden^um 
in the mass- range A=90-106, Whereas the spectra in A<96 
have been known for quite some time now, the spectra for 
A>98 nuclei were investigated relatively recently by Cheifetz 
et , al^~ 

90*- 10 2 

A striking feature of the spectra in “ Zr is the 

sudden onset of deformation at N = 60; whereas one observes 

large energy separations (>_0 # 92 MeV) between the ground and 

90—98 

the first excited states in Zr, the yrast spectra (with 

= 3 + ) in are almost rotational with 1 i 

max i 

i 

! 

I 



4 


* 
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92 94 96 98 100 102 


) Experimental spectra of the even Zr isotopes. 




94 96 98 100 102 104 106 

Kb) Experimental spectra of the even Mo isotopes. 



CE^ + - E q+ ) = o,21 MeV. In the doubly-even Mo isotopes, on 
the other hand, one observes a distinct but gradual increase 
of deformation in going from A = 92 to A = 106, 

The mass region A = 90-106 offers a nice example of 
shape transition. The region includes at the one end nuclei 
which are easily describable in terms of sirrple shell-model 
configurations?' 3 At the other end of this region (A>100) 
one comes across unambiguous evidence of rotational collecti- 
vity. 

Microscopic calculations of the yrast spectra in the 

Zirconium region have not been atterrpted so far. However, 

some attempts have been made in the last few years with a 

view to identify some possible mechanisms that could be 

held responsible for the dramatic onset of large deformations 

4 7 

at A = 100, In this connection, Federman and Pittel ' 
have suggested that the onset of rotational features in 
the highly neutron-rich Zr isotopes may be crucially linked 
with the deformation of the Z = 40 proton core resulting 
from ^■l < 3g/2^ “ ^ % /2^] T=0 interactions. Their calcula- 
tion of the deformation energies 5 - and these were corrputed 
only for the Mo isotopes — also gave some indications 
regarding a possible connection between the onset of deforma- 
tions at A = 102 and the T = 0 interactions in the 

(ig 9/2 - l9 7 /2 ) 803 (lff 9/2 ' lh ll/ 2 > • 

94 

The above calculation was performed assuming a Sr 
inert core and errploying as an effective interaction the 
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surface delta interaction. The active neutrons and protons 
were restricted to only the following orbits outside the ^Sr 
core; 2P,jy 2 • 1 "%/2 an< ^ ^5/2 for P rotons and ^ s i/2 # ^ 3 / 2 ' 
197/2 an< ^ l^n /2 f° r neutrons. 


Employing the above valence spaces for protons as well 

as neutrons, Pederman, pittel and Campos^ have also carried 

98 

out shell model calculations for the nucleus Zr with a view 

100 102 

to correlate the onset of deformation in * Zr with the 
increased involvement of (1 Qn/o^ 1 configurations. 


/2 


In their shell model calculations the first-excited 0 level 
98 

in Zr (with E x = 0.9 MeV) turns out to possess a wavefunction 
spread over many configurations. It is seen that the extent 
of configuration mixing in the 0 ^ state depends very much 
on the relative single-particle energies of the 1 ^ 9/2 and 
ig 7 ^ 2 orbits. Having established a connection between the 
simultaneous occupation of the (lgg^)^ and (lg 7 / 2 ) n orbits 
and the "collective" nature of the first-excited state in 
^Zr, Pederman et al f have conjectured that the rotational 
states in 100 Zr may also be possessing a structure quite 
similar to that of this state. 


An iirportant feature of the above calculations is the 
94 

assumption that Sr is a good core. The purpose of the 
present work is to examine the deformation systematics in 
the A = 90-106 mass region by employing a valence space large 
enough to permit an unbiased study of the shape transition. 

To this end we employ pairing-plus-quadrupole-quadrupole 



interaction operative in a valence space spanned by the 

2p l/2 ' lg 9/2 ' 2d 5/2 ' 3 s 1/2 ' 2d 3/2 ' lg 7/2 and lh ll/2 orbits 
for protons as well as neutrons. The nucleus ^Sr has been 

considered to be an inert core. 


The present calculation in the extended space is quite 

successful in reproducing the observed yrast spectra. A number 

of interesting features associated with the dramatic shape 

transitions in the Zirconium region are revealed. The results 

indicate unambiguously that the observed onset of large 

deformations in this region is not compatible with the 
4 , 04 

assumption of an inert Sr core. Further, it turns out that, 

5 

in contrast to the suggestions made by Federman and Pittel , 
the observed rotational spectra can be quite successfully 
explained without invoking strong T = 0 interactions for the 
( lgg/2 "* lg 7/2^ configuration. The onset of deformation and 
the coexistence of deformed and spherical shapes in the 
A = 90-106 region are seen to be related mainly to the 
involvement of the lh.^^ level, which possesses a large 
quadrupole content, in the valence space. 


In section III. 2 we discuss the choices of input para- 
meters, such as the spherical SPE's and the pairing-plus- 
quadrupole-quadrupole effective interaction that we have 
employed in the present work. Section in, 3 describes the 
Cranked Hartree-Fock (CHF ) approximation which has been 
employed in carrying out some preliminary studies of deforma- 
tion trends in the Zr isotopes. The CHF results in the space 



(3s-2d-lg-lh 11 ^ 2 ) are discussed in section III. 4. In section 
III. 5 we discuss the Cranked Hartree-Fock-Bogoliubov (CHFB) 
prescription for studying the yrast spectra. In order to 
examine explicitly the role of protons vis-a-vis the onset of 
deformation in Zr isotopes, we have carried out CHFB calcula- 
tions employing valence spaces that include also the 2p.jy 2 
orbit in addition to the orbits included in the calculations 
discussed in section III. 4. The results are presented in 
section III. 6. Finally section III. 7 contains some concluding 
remarks. 


I II . 2 The Input Parameters of the Calculation 

The spherical SPE*s that we have used in our 
calculations are (in MeV) : e(lg g ^ 2 ) = 0.0, e (2d 5 ^ 2 ) = 5.4, 
f( 3s 1 / 2 ) = 6.4, e (2d 3/ , 2 ) = 7.9, e (lg ? ^ 2 ) = 8.4 and = 9 . 0 . 

This set of input SPE's is 'exactly the same as the one employed 
in a number of successful shell— model calculations by Vergados 

P 5 

and Kuo as well as by Federman and Pittel except for a slight 
lowering of the ener< 3Y b Y °- 4 MeV - 


The two-body effective interaction that we have employed 
is of "pairing + quadrupole— qua dr up ole" type. The quadrupole— 
quadrupole part of this schematic interaction has a tendency 
to give rise to the quadrupole deformations in nuclei and the 
pairing term favours sphericity since it acts only between a 
pgf £• of nucleons in time— reversed orbits. The pairing part 



(III. 2.1) 



Here 'a' stands for the set of quantum numbers (nl jm) . The 
state a is the same as a but with the sign of m reversed: a 
= (nlj-m). 

Further, 


s a = (-1) 


j-m 


(III. 2. 2) 

The quadrupole-quadrupole part of the two-body interac- 
tion is given by 

+2 

v a CT = ~ -2 X I I lY><5|q 2 u | B >(-l) V cjcjc C (III. 2. 3) 

4 a3y<5 y=-2 “ ^ c Y 


2 

The operator q^ is given by 


q 2 = UST/S^rV (e,jzf) 


(III. 2.4) 


The pairing-plus-quadrupole-quadrupole interaction 
leads to absurd consequences if it is allowed to operate in 
a complete set of states. Consider the pairing part of the 
interaction in the coordinate representation: 


< r‘+ / r*+|V p | r+, r + > = - G 


(III. 2. 5) 


The particles must be initially both at r and finally both 
at r 1 since they must be in time-reversed states. The 
constancy of the matrix element (III. 2. 5) - it has the value 
-G for arbitrarily large |r-r' ] - is clearly quite unphysical. 
The usage of this interaction is, however, justified for 
restricted spaces where it can simulate the effects due to 
the short-range components of the interaction. The quadrupole. 
quadrupole part also yields, via a Hartree type calculation, 
an effective quadrupole operator which is unbounded in a 



complete set of oscillator states. Here again the interaction 

derives its validity from the limited size of the space. The 

optimum size of the oscillator space is usually two major 
9 

shells • Kumar and Baranger have shown that the matrix 
elements of the q. q force computed for two major shells are 
quite similar 'to the ones obtained from Hamada- Johnstone • 
potential, using the G-matrix plus core polarization. 


III. 3 The Cranked Hartree-Fock (CHF) Method 

In this section we present a brief outline of the 

CHF prescription that we have employed in the next section 

(subsection III. 4. 3) in order to examine some questions concern- 

ing the coexistence of the spherical and deformed states in 

l°0,102 zr _ 


10 

It was shown by Beck, Hang and Ring and Villars and 
1 1 

Schmeing-Rogerson that the energy corresponding to a parti- 
cular angular momentum (J) state projected from an intrinsic 
state 0 which contains predominantly a K=0 (K=<£f \j^0 > = the 
projection of J along the body-fixed Z axes) state and small 
K =£ o components is approximately given by the following 


Ej 5 - 


% I > 


«HhU> - < <f> 1 Jy ! 4>>/l 


+ J+i) ~ < ~ 


(III. 3.1) 


where 


(Jj* S 


< HJ x > " <H><J X > 

< > 
y 


5 ' 


(III. 3.2) 


/ 



Here the brackets < > imply that we are taking the expectation 
values with respect to the intrinsic state $ . 

Now, carrying out HP type variations on Ej is equivalent 
to the cranked HP problem 

6 [ <H — a> J > J = 0 (III. 3.3) 

X. 

with the subsidiary condition 

<J X > =4j(j+l) (III. 3.4) 

provided the following two conditions are satisfied. 

2 

(I) The term which can be considered as a zero-point 

contribution to rotational energy, has no significant J 
dependence and hence can be treated as a constant during the 
variation on Ej ; and 

(ii) The frequency cu appearing as the Lagrange multiplier 
in eq. (III. 3.3), which is needed to satisfy the subsidiary 
condition (ill. 3,4), is the same as the to* of eq. (III. 3.1). 

The condition (i) has been tested and found to be valid 

10 

for some values of « by Beck, Mang and Ring . The condition 

12 

(ii) has also been tested by Sharma et al. in some 2s-ld 
and 2p-lf shell nuclei and found to be approximately valid. 

In the present work we have solved the eq. (III. 3.3) 
with the constraint (III. 3. 4) for some Zr isotopes by employing 
pairing-plus-quadrupole-quadrupole effective interactions in 
the 3s~2d-lg-lh^\jy2 


space, 



III « 4 * 1 A Preliminary Study of the Deformation Trends 
in the Doubly-even zirconium and Molybdenum 
Isotopes 

We discuss here first the results of our deformed 
HF calculations carried out by employing the pairing + q. q 
interaction. The strength • x (appearing in eq. (ill. 2. 3) for 
like-particles was taken as x nn = -0.02 MeV b“ 4 for the Zr 
isotopes. Here b is the oscillator parameter. The intrinsic 
states in Mo isotopes were computed by setting Xnn = )( = 

-0.016 MoV b~ 4 and X np = -0.0352 MeV b“ 4 . In all these 
calculations the strength of the pairing interaction was fixed 
(through the approximate empirical relation G = 18/A) at 
G = 0.18 MeV. 


The present calculations in the HP framework, in which 
the pairing effects have been ignored, are intended' to serve 
only as a basis f or a qualitative discussion of the deformation 
systematic s in the Zirconium region. We shall see later, in 
section III. 6, that the results are not drastically changed 
when we include the pairing effects. 


90, 96 

Wo first discuss the nuclei Zr. The observed spectra 

in those nuclei are characterized by large values (v2 MeV) 
of (E - E ) sepa-t 3 ^ 0113 * This can be interpreted in terms, 
of sub-shell closures for the lg ^2 an< ^ ^d- 5/2 orbits. The 
present calculation reproduces this feature - the minimum- 
energy HE' solutions that we obtain for 90/96 Zr with the present 
set of input parameters are indeed spherical with the structures 
!lg, /2 )“ and [(13 9 /2 )10 (2d 5/2 )S ]j=0 respectively. 
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The present approach differs significantly from that of 
earlier scuuies in that a light core (the nucleus ^^Zr) is 
assumed* However it is essential when using such a light core 
to first confirm that the calculations do reproduce the experi- 
mental 11 single-particle" spectra in the ^Zr-plus— a nucleon 
system* In Table 111*1 we have given the ^Zr SPE f s resulting 
from the relevant spherical HE solution* The calculated SPE 1 s 
for various orbits compare favourably with the values obtained 
from the (d,p) as well as (p,p') reactions . 


TABLE III.-l 
91 

Comparison of the calculated Zr SPE's with their 
experimental values 


Neutron orbital 

Energy 

Theory 

(MeV) 

Expt* 

2d 5/2 

0.00 

0.00 

n 1/ 2 

1.04 

1.20 

2d 3/2 

2.30 

2.08 

** 7/2 

2.90 

2.12 


92 94 98—102 

We next present the results for the nuclei ' ' Zr 

and 92 - 106 Mo in Table III. 2. Our calculation yields small 

<0 2> values for the nuclei 92 ' 94,98 Zr. This agrees with the 
o 

observed lack of rotational collectivity in these nuclei. 

2 

Note the sudden jump in the intrinsic < Q q > values m going 
from 9 ®Zr to 299 Zr; whereas the < Q 2 > value in 98 Zr is 18.74 b 
out of a maximum value of 62.7 3 b 2 for the given space, that 
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TABLE III. 2 

The calculated values of the intrinsic quadrupole moments (<Q 2 >) 
and the sub- shell occupation numbers for the Zr and Mb isotopes. 
Here < Q 0 > mox l^ves the maximum value of intrinsic quadrupole 
moment for each isotope, The quadrupole moments have been 
computed in units of b 2 . 


Isotope 


< q 2 > 
w o 


< Q 2 > 
o max 


Sub- shell occupancy 


14.5 


17.0 


51.7 

55.6 


9.94 

9.96 


1.49 

3.25 


0.01 


0.06 


18.8 


62.7 


9.99 


5.53 


0.65 


32.0 


62.6 


9.76 


3.89 


1.82 


38.9 


60.6 


9.65 


3.59 


2.62 


12.7 


63.6 


9.84 


0.14 


31.9 


67.6 


9.66 


1.13 


0.14 


41.4 


71.6 


9.61 


1.66 


1.24 


45.7 


75.6 


9.78 


2.57 


1.60 


53.6 


78.7 


9.73 


2.41 


1.76 


59.8 


78,6 


9.68 


2.31 


1.87 


59.5 


76,6 


9.76 


3.27 


2.82 


62.9 


74.6 


9.74 


3.26 


2.84 
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for th>_ nucleus Zr is 31.97 b^ which is roughly 50 percent 
2 

of the < Q q > max in this nucleus. The present results, therefore, 
explain, at least qualitatively, the experimentally observed 
dramatic onset of large deformations at "^^Zr. 


The results for the nuclei ^ ^^Mo also reproduce 
successfully the deformation trends that are implied by the 
observed spectra. The nuclei ^^Mo display well— developed 
rotational spectra. This is easily explained qualitatively 
in the present calculation where we obtain large intrinsic 
deformations - roughly 80 percent of the maximum possible 
values - for these nuclei. 


The reason for the sudden increase in intrinsic deforma- 


tions at N = 60 now becomes immediately obvious if one considers 

the subshell occupancies for the 9i v e n in column 7, 

100 102 

Table 13 T.2. The nuclei Zr and Zr have two and four 


neutrons in the orbit lh respectively. Whereas the sudden 

98 100 

change in the Ih^^ occupancy in going from Zr to Zr is 

seen to lead to an increase of about 65 percent in the computed 
2 

< Q > values, a further change of the occupancy in going from 
100 102 

Zr to Zr results in an increase of about 20 percent. 

The effective Nilsson deformations for the nuclei 100 'i° 2 Zr 


consistent with the computed lh^^/2 occu P anc ^- es are shown 
in figure 2. 


The results for the lhjj/j sub-shell occupancies in the 
Mo isotopes further emphasize the role of the lh ^^2 orbit 
in producing large deformations — and the associated 



Comparison of the deformations of the HF states of the 
N = 60 and 62 isotopes of Zr. The values of k = <j 2 > are 
given on the right.. 
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quasi rotational features - in the A = 90-106 region. In view 
of the fact that Ho already has six neutrons in the 

1/2 otbit, one expects a decrease in the <Q^>* values in 
the heavier Mo isotopes. 


The calculated results for the sub— shell occupancies for 

the l9g^2 / 2c %/2 anc ^ ^7/2 or kits in the Zr and Mo isotopes 

(see columns 4-6, Table III. 2) indicate that the nucleus 98 Zr 

can approximately be considered an inert core even in the 

quasi rotational isotopes 10 °- lo2 Zr and 102 ~ 106 m O ; the maximum 

depletion of the orbit is seen to be only about 3 percent 

in these isotopes. The 2d 5/ , 2 occupancies in 98 “ 106 Mo are, 

however, seen to range only between 2.57 and 3.26. This shows 

4-6 94 

that the usual assumption of an inert Sr core is not 
valid for nuclei with N>56. 


1 1 1 , 4 . 2 The Role of T=0 Interactions in Producing 

Deformations in Molybdenum Isotopes 

5 7 

In the earlier work ' in which the surface 

. 94 

delta interaction was employed in conjunction with Sr core, 
the T = 0 interaction had turned out to be the major source 
of deformations in the Mo isotopes. In order to assess the 
role of this component of the interaction in producing deforma- 
tions in the context of present calculation which involves an 
extended space, we have obtained the HF states for 100 “ 106 Mo 
with the full interaction as well as the T = 0 part alone. 

The results are given in Table III. 3. 
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TABLE III. 3 

2 

The < Q q > values in some Mo isotopes calculated with the full 
interaction ( V j= ul ^ as well as its T = 0 component ( v T _o) 


V full 

V T=0 


100 MO 

102 M 

Mo 

53.59 

59.79 

09.78 

18.59 


104 MO 

o 

O 

59.49 

62.86 

22.58 

21.00 


We find that the T = 0 part plays only a minor role in 
producing deformations in the Mo isotopes. This relative 
suppression in the deformation producing role of the T = 0 
component in the present case is closely related to the fact 
that the choice of 80 Zr as a core leads to a reduction in the 
numbers of n-p pairs in corrparison to the n-n and p-p ones. 


111.4,3 The First-excited 0 + Levels in th e Doubly- even 
Isotopes of Zirconium 

A significant feature which characterizes the 
experimental spectre in the nuclei 90 - 100 zr is the occurrence 
o£ an excited 0 + level at rather low energies. In Zr the 
first 0 + level occurs at 1.59 MeV, at roughly the same energy 
as in the nuclei 90 ' 94 Zr. In 98 Zr, however, this level appears 
at 0.85 MeV and in 10 °Zr it suddenly drops to just 0.33 MeV. 

Federman, Pittel and Carrpos 6 have recently tried to 
examine the nature of the first-exoited 0 + level in the 
nucleus 98 Zr in the framework of shell-model by considering 



, i 94 

the nucleus Sr to h<=> 

an inert core. Their calculation 
indicates a "collective" nature for. this state in the sense 
that it is made up of a large number of shell-model configura- 
tions. However no calculations have so far been carried out 
to study, in a systematic manner, the (E y - E +) trends for 
the nuclei zr. 

We present here an HF-oriented study of the excitation 

energies of the first o' 1- levels. In Table III. 4 we give the 

results concerning the position of the first excited 0 + levels 
9 2 1 0 2 » 

in the nuclei Zr. The 0 + positions in the nonrotational 

. 92-98.. 

nuclei Zv have been obtained by corrrputing the difference 

in energies between the minimum— energy HP state and the lowest 
1 particle- 1 hole state with K = o * The excited 0* levels 
in the qua si rotational nuclei 100 ' 102 2r have b een calculated 
Ly obtaining the non— axial/ excited HF solutions* It may be 
mentioned here that an explicit inclusion of the pairing 
correlations is not likely to affect our estimates for the 
(0 -0 ) separation. This is because the pairing correlations 

are expected to lead to a lowering of the two HF intrinsic 
minima (or a HF state and a Ip-lh state) by approximately 
equal amounts. 


TABLE III, 4 

Corrparison of the calculated as well as the observed positions 
(in MeV) of the first excited 0 + states in the Zirconium isotopes 


92 

Zr 

2.27 

1,40 


94 

Zr 

2.22 

1.32 


96 ^ 

Zr 

2.14 

1.62 


98 „ 
Zr 

1 . 66 

0 f . 86 


100 zr 

0.53 

0.33 


102 Zr 

1.62 

'1.50 




H At 






Calculated 

Experimental 


The results concerning the 0 + ' positions are in good 

qualitative agreement with the experiments - the sudden drop 

in the excited 0 position in ^®Zr, after remaining virtually 

constant in Zr, is reproduced. In 100 Zr the excited 0 + 

l.ie.j at 0.53 MeV in the calculated spectrum, in reasonable 

agreement with the experiments. This o"^* state is seen to arise 

from a weakly- deformed intrinsic state (with <q 2 > = 16.8 b 2 ) . 

In view of the highly deformed nature of the ground state band 
100 _ .. . 

xn Zr, thxs nucleus appears to provide an example of the 

coexistence of the shapes with widely varying degrees of 

deformation. As discussed earlier, the addition of two more 

1 ^ 11/2 neutrons to ^"Zr leads to a significant increase in 

the intrinsic deformation of the yrast levels. This leads 

to a further lowering of the yrast 0 + state relative to the 

excited one since the latter has a nearly- spherical parentage* 

Thus the excited 0 + lies at 1.62 MeV in the calculated spectrum. 

Importantly enough, no excited 0 + has so far been reported 
10 ? 

in Zr at least upto 1.55 MeV of excitation. 

In figure 3 we display the theoretical yrast spectra 

100 102 

that we have conputed for the deformed nuclei ' Zr in 
the framework of the CHF procedure discussed in section III. 3. 
Various non-axial intrinsic states, <i>(w), were obtained by 
solving the equation III. 3. 3 — »and this is just the HF 
equation for the Hamiltonian (H-a»J x ) —for various values Of to. 
The intrinsic states for various yrast levels were then 
obtained by subjecting various $0*0 1 s to the constraint III. 3.4. 
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yrast spectra 


the calculated and the exp*r imen,al 
in nuclei ,0S Zr and ,0 *Zr. 
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The results for the nuclei iOO/lC^ 

are presented in 

Table III. 5. . 

It is interesting to see th »+- fK a 

y u oee that the present CHP calcula- 

tion, which is rather elementary in the sense that it ignores 

pairing effects as well as the participation of the 2=40 core 

protons, leads to the comnuteri . 

“pu ed spectra which are m reasonable 

ngreoment with the experiments. The calculation reproduces the 

observed positions of all the yrast levels, except the j* =6 + 

level in 2r, within an accuracy of 140 kev. 

1 r X . j , 1 Jhe Cranked Hartree-P o ck-Boaoliubov ( CHFBi Method 

In this section we discuss the general CHFB 

I .. ; ,.ot. / which i.i applicable for any given two— body effective 
14 

m to ruction . 


Con cider the many-body Hamiltonian 
H » l <ileli>a; a. -i- 4 E <ij|v |kl> a T at a,a, (III. 5.1) 


+ ^ + 7 S -i hr I VI ^ ~ + - + 

1 4 ijkl 


£, < ije|i>a i a i + j l . ,. V1 '--j i - a i — - 


in iicoc "ijkl" span the active valence single-particle 
nt-iLes contained in the model space, and at and are the 
particle creation and annihilation operators respectively. 
The effective nucleon-nucleon interaction V appropriate for 
the : no del space is anti symmetrized, so that 


<ijjv a |kl> = <ij |v|kl> - <ij [V|lk> (III. 5. 2) 

The main idea involved in the CHFB or HFB theory is to 
obtain a transformation from particle -co-ordinates to qua si- 
pur ti cl o co-ordinates such that the quasiparticles are 
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relatively weakly interacting. That is we want to express H 
as H = w -L. U . u (III. 5. 3) 


H = E + H + H 

° <2P qp-rnt 


where E q is the energy of the quasiparticle vacuum, H 

describes the elementary quasiparticle excitations, and 

^qp~int c h arac "k er izes what we hope is a weak interaction 

between quasiparticles, in HFB this last term is neglected, 

so that H is approximated by an independent quasiparticle 

Hamiltonian. It is obvious that the neglect of H . will 

J qp-rnt 

result in a loss of symmetry - the Hamiltonian (H + H ) will 

o qp 

not possess, in general, the same symmetries as the initial 
Hamiltonian H. 


Just as in the BCS theory, HFB includes pairing, correla- 
tions by introducing number nonconserving wave functions. The 
Hamiltonian H is replaced by 


H* = H - X N - X N 
p p n n 


(III. 5. 4) 


where the Lagrange multipliers X^ and X n turn out to have 
the physical interpretation of proton and neutron Fermi 
energies. The Lagrange multipliers are adjusted so that the 
number operators 

N p = j[ a ip a ip ' N n = J a in a in (HI. 5. 5) 


have the expectation values 

<$ o lN p 1$ o >= Z ' < ^o^n^ $ o > = (III. 5.6) 

where |$ Q > is the quasiparticle vacuum. 



The independent quasiparticle Hamiltonian is, in general, 
not rotational ly invariant. Therefore angular momentum is not 
a good quantum monger for j $ q >. The Inglis cranking model 15 " 16 
is invoked to approximate angular momentum projection. 

Define 

H ' = H - x p N p - *n K n - “ J x (III. 5. 7) 

where the Lagrange multiplier, e is adjusted so that the 
angular momentum component 

JL 

J x = } _ <i | J x |j > at aj (III. 5.8) 

has the expectation value 

<* 0 |j x ! V = (III. 5.9) 

The expressions given above correspond to the intuitive 

picture of a nucleus that rotates with angular velocity a 

17 

about the x axis. Mang has shown that the cranking prescri- 
ption also results when one considers an appropriate approximate 
treatment of angular momentum projection. 


Although the mean value of angular momentum in the y and z 
directions vanishes. 


<» 0 l J yl V = < ' t ol J Z f»0 > = °' 


(III. 5. 10) 


there may be significant angular momentum fluctuations in 
these directions. The cranking term (-aJ x ) violates time 
reversal, which technically complicates the calculation of 
the HFB intrinsic state. 
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111.5*2 Derivation of the CHFB 


lation* 


We next demonstrate how the schematic separation 
of the Hamiltonian/ indicated in equation (ill. 5. 3), can be 
realised. A quasiparticle is an excitation of the system. In 
the vacuum state appropriate to HF theory/ all single-particle 
levels | a> below the Fermi energy e are occupied, and all 
levels above e j, are empty. Thus in the HP theory the quasi— 
particle creation operators q* are defined by 


U a > 'P 

^ £ a < 


(III. 5. 11a) 
(III. 5. 11b) 


where (111.5.11a) signifies creation, of a particle and 
(111.5.11b) represents creation of a hole. 


In BCS theory the quasiparticles are defined by the 
special Dogoliubov transformation 

<£ “ a t * V a a !» (HI. 5. 12) 

where loT> is the time reverse of j a>, and U a and V a are 
real. 

The HFB theory employs the general Bogoliubov transforma- 
tion q+ = Z (U ±j a+ + V ±j aj) (III. 5. 13) 

in which each quasiparticle is a linear combination of all 
particle creation and destruction operators. If there are N 
single-particle basis states, then U and V are N x N complex 
matrices, and q and a are N- component vectors. 

The 2N x 2N linear transformation 



(XIX. 5. 14) 


Ki 

| 

U V j 


r 4.1 

j 

1 = ; 

is JL. 

, v u ■ 


a i 



r 



is required to be unitary. Equations (III. 5. 14) can then be 
inverted to express the particle operators in terms of the 
quasiparticle operators 

a i = j (U ji <4 + V.. qj ) (III. 5. 15) 

The unitarity constraints are 

UU + + W + = U + U + VV* = I (III. 5. 16) 

UV * VU = U + V + W* = 0 (III. 5. 17) 

.t. 

where U denotes complex conjugate and U denotes the trans- 
pose of U. These equations may also be derived by requiring 
the particle and quasiparticle operators to satisfy the fermion 
anti commutation relations 


[®i ' a jj“ *ij ' 

[ a i * a j] =f a i' a j 1= 0 

(III. 5. 17) 

5 ij 

- <3j] = 0 

(III. 5. 18) 


The quasiparticle vacuum ]$ o > is defined by the requirement 

q. |* o > = o (for all i) (III. 5. 19) 

Prom eq. (III. 5. 18) it follows that the solution of eq. - 
(III. 5. 19) is 

|$, >=N n q i |0> (III. 5. 20) 

° i 

where N is a normalization constant and (o> is the particle 


vacuum 





The density matrix p and the pairing tensor t are 
defined by 

p ij = <4 ol a j 

t ±i = <» 0 |a j a 1 |j 0 > 

It follows that p is hermitian and t is antisymmetric. 

Also t is nonzero only for number nonconserving wave functions. 
To evaluate p and t, substitute eg. (III. 5. 15) into eqs. 

(III. 5. 21), (III. 5. 22) and apply eqs. (III. 5. 17), (III. 5. 18). 


(111. 5. 21) 

(111. 5. 22) 


The result is 


P = V + v 
t = V' : " U 

The constrained Hamiltonian is defined by 

H' = H - C 

c = x p N p + x n N n + “ J x 


(111. 5. 23) 

(111. 5. 24) 

(111. 5. 25) 

(111. 5. 26) 


where H is given by eq. (III.5.1). By substituting eqs. 

(III. 5. 21), (III. 5, 22) into eq. (III. 5. 25), _H can be expressed 

as 


H ' = (III. 5. 27) 

where liL contains n uncontracted operators, and 

(III. 5. 28) 

' 2 ' " ' iJ 

I I 

+ 


n 


H* s Ir[(e - C + r -*- ) P At + 1 

* *Q W ' 4t- 


4 - £ (X - c) i 3 ! 4 a r + $ ij 4 a r 

£ 4j : a i a j ! 


(III. 5. 29) 



bl 


= * i/kl <1J l V Al kl>: 4 ai V (III. 5. 30) 

Ihe HP Hamiltonian X f the HF potential r and the pair potential 
A are defined by 


r 


ij 


= e + r 
= k \ 


A 

ij 



I <ij 
kl 


l V al kl>t kl 


(111. 5. 31) 

(111. 5. 32) 

(III. 5.. 33) 


It follows that and p are hermitian and that A is anti- 
symmetric. For any normal product :0:, the expectation 
value < $ o ) :0: | ® o > vanishes by construction. Consequently 

, » i 

< $ 0 1 H 2 + l $ o > = 0/ so ’ t ^ at tlle vacuum expectation value 

i i 

< I H ] $ Q > is H . The quasiparticle vacuum energy is 


<$ C |H| *0 > " <<b o\ H + c I ® 0 >= Tr[( « + | r ) P+-|At + ] (III. 5. 34) 

The first term on the right-hand side is the HF binding energy, 
and the second term is the pairing energy. 


The ^ua sip article excitation energies are obtained from 

i i 

Hg. The remainder H^ describes the quasiparticle interactions, 
which are neglected in HFB. 


I 

As mentioned earlier our aim' here is to express H as 


H = Ho + I E. q. q. + 


(III. 5. 35) 


This will be achieved if acquires the form of an independent 
quasiparticle Hamiltonian 



? E i%% 

i 


(III. 5. 36) 



[ H 2 ' 4 ] ■ 


we 


Inserting' eq. (111,5,36) into the commutator 
have 

[ H 2 , q+ J = E qt = E. z (U a"!) + V. . a .) (III. 5. 37) 

+ ± X JJ X JJ 

Alternatively, using eq. (III. 5. 29) we obtain 


[ H 2 - qp = I t < X - C) 1k U ik + A 1k V. k ] a + 




jk ik T *'jk ik j 


+ I [-slu. v - (X-otv.Ja, (III. 5.38) 


jk 


jk ik 


jk xk J ~j 


Equating the coefficients of at and aj in eqs. (III. 5. 37) and 
(III. 5. 38) we have 


fx 


c 

* 

A 


(X. C) 


* 

f tn 1 


U 1 


. X . 


U. 

j ' 

1 ->■ 

V. 

1 a-j 

= E. 
x , 

l 

^ij 


(III. 5. 39) 


where U. denotes the vector (Uh ^ , U j_2*** U ilP and N is the 
number of single-particle states and similarly for ^ . These 
are the HFB equations. 

We now discuss several simplifications which result 
provided one exploits certain symmetries associated with the 
J x operator. We also make use of the simplified nature of 
the interaction. v • 


III. 5. 3 Symmetries 

It is convenient to separate the harmonic oscilla- 
tor states into two sets. The first set contains the states 
jk > , which are now restricted to have (m - ■§) equal to an 
even integer. The second set contains the time-reversed 
states jk > * T jk> , which have (m - ■§) equal to an odd integer. 



The phase convention is T|nljm x > = C-l) j - m+1 1 nl j-m T ►. For 
quasiparticle operators of the form 


q'. h 

X 


+ 


l (u ik < + v lh a T?) 


l <U ik + v ik a k ) 


(III. 5. 40) 

(III. 5. 41) 


P does not connect the set of states |k> and |k>. For t, 
however, only the matrix elements connecting |k> and |k > are 
non— zero. Since the nucleon-nucleon interaction conserves 
angular momentum, the average potentials (r,A) follow the 
partitioning of the densities (P,t). That is 


p = 

r \ 0 5 

t = 

f° h; 


'o 0 p 2 > 

/ j 

ft 

to 

O 

V- 


(III. 5. 42) 




(X ■ 

1 

0 1 

1 / 


0 


X 

as* 




A = 





K 0 

V 



^9 

0/ 


where n and X are hermitian. Substituting (ill. 5. 42) into 
the HFB equations demonstrates that the energy matrix decorrpo- 
ses into two blocks and that the eigen-vectors retain the form 
of esqs (III. 5.40), (III. 5. 41). Consequently the partitioning 
propagates through each iteration. 

Time reversal symmetry is irrposed by requiring that 


4 = t 4 T ' x 


no that 


^ik 


* V., 

u ih ’ lk 


* 

“ V ik 


(III. 5. 43) 
(III. 5.44) 
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Then the quasiparticle vacuum is time-reversal invariant, and 


P-3 = p" .X _ X t + - t A + _ A 

2 y l * 2 1 * r i “ t i ' A- = A 1 


(III. 5. 45) 


Consequently onlv the block of the energy matrix corresponding 
to oq* (III# 5 # 40) need be diagonalized 


■,** V ✓ 

"l : hi' 5. 

(\ - c) , tv. I ‘ E i v 1 
1 / 1 1 


(III. 5. 46) 


Prom eqs (III. 5. 45) it follows that the time reversed 


quasip articles q* and qr are degenerate 


E I = E i 


(III. 5.47) 


Oy using eq. (III. 5. 45) along with the unitarity 
constraints on t^ and , we find that and t^ commute, 
Go a basis }aa> exists, where 


a a ” £ D ak a k 

- * 4 4 


(III. 5. 48) 


(III . 5. 47) 


;;uch that p and t have the canonical form. The constraint 


2 + 

p - p = tt implies that 


:[p„ (1 - P„)l ** 


(III. 5. 50) 


which further indicates that the HPB quasiparticle vacuum 


may bo expressed in the BCS form 


I p > = n (U + V a + afi 1 0 > 

* v O a a a a 

a 

The paired orbitals are now related by time reversal- 


(III. 5.51) 
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A wave function is defined to have triaxial symmetry 
if it is invariant, except for a multiplicative phase factor, 
with respect to rotation by ir about each of the three principal 
axes and rf lection through each of the three principal 

planes ce , where 

R^(tt) = exp (-iffCR) (111.5,52) 

°i = p R j_ ( ir ) (III. 5, 53) 

More i denotes x,y, z and a is reflection through the yz-plane. 

A 


The effect of R^(ir) on the single-particle state jjm>is 
given by 

R v (*)|jm> = e - iffj |j-m> (III. 5. 54 

Ry (if) |jm> = (-1) |j-m > (111-5,55 

R (ir) |jm> = e“ i7rm |jm> (III. 5. 56 


R (it) |jm> = e” 1 "™ 


(III. 5.54) 
(111-5*55) 
(III. 5. 56) 


VA ncc | k > is restricted to states with(m - ^) even, the 
particle operators transform under R z (ir) as 


\ Kt R z 1(7r) = - 1 


(III. 5,57a) 


R„ («) a ^;; R z 1 ^ ir ^ = +i 


(III. 5. 57b) 


Consequently the quasiparticle operators q+ and qj are 
invariant with respect to R^U) end therefore so is the 
vacuum |« 0 >. Also notice that a y | jm > = T |jm>. so that 
and T are equivalent if the quasiparticle coefficients V and 

V are chosen to be real. 
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lo oummari ze, if the qua sip article operators are chosen 
as 

q± ” k <U ik + V ik (III. 5.58) 

^ " l (U ik ^ - V 1X V (III. 5. 59) 

where U and V are real and (m^. - —) is even, and. if each quasi- 
particle operator contains particle operators possessing the 
same parity, then the HFB vacuum |$ 0 > is time- reversal invariant 
and tri axially symmetric. 

The self-consistent symmetries discussed above are the 

i 

ones satisfied by H. The constrained Hamiltonian H 

(uq. III. 5. 4) violates most of the symmetries mentioned here* 

For instance J anti commutes with the following operators: 

*A> 

time reversal, rotations about y and z axes and reflection 
through the xy and xz planes. Therefore cranked HFB wave 
functions cannot display time reversal, triaxial or axial 
symmetry. This loss of symmetry makes the calculations very 
difficult. 

The symmetries preserved by J x are parity, reflection 
through yz plane, and rotation of tt about the x— axis. Since 
°x ~ PR x iv) ' only tWO ° f these symmetries are independent. If 
the conventional | jm> basis is used, then the reflection 
symmetry a does not reduce the dimension of the CHFB equations. 
This is because J v is a' linear combination of raising and 
lowering operators, so that J x is not block diagonal in the 

quantum number m* 



To take advantage of the symmetry cr x , Goodman introduced 
a single particle basis that simultaneously block diagonalizes 
the J x term and the HP and pair potentials. Since a 2 = ~1 

X 

when acting on on^—f ermion states* the operator a has only 
two eigenvalues, -i and +i, with corresponding sets of eigen- 
vectors that we denote by ]k> and |k>. Because a and J are 
commuting normal operators* it follows that there are no non- 
zero matrix elements of J between eigenstates of cr belonging 
to different eigen values of a . That is, 

<K| J X |K > = 0 (III. 5.60) 


The eigenvectors are easily determined, since a in the two- 

X 

dimensional basis |k k >has the representation 

o x = -i (° J) (III. 5.61) 

so that the eigenvectors are 

| K > = 2“^ [ J k > + fk > ] 

(III. 5. 62) 

I k > = a -* 5 [_| k > + |k>] 

Notice that | K > = t|k> . Since J x and T anti commute it follows 
that 

<k|jJk'> = -<k|j!k‘> (III. 5. 63) 

* X "X. 


nqs. (III. 5,60) and (III. 5.63) indicate that J x is block 
diagonal in the cr x basis Jk 1 K 2 ... K N K i K 2 ••• ^ > : 


'x 



y 


xj 


0 

-j. 


(Iil.5.64) 



If the quasiparticle operators are chosen as 
^ = R (D iK 4 + V iK ^K> 

3 - k 4 + ^iK -K> 

it can then be demonstrated that the CHFB equations reduce to 
the form 

d - “ j x \ 

- (Xj + »J X ) 


* + 
l 


\ 

f-v } 


U. 


l , 


-»■ 


v. . 

) 

V 


' % 

u 


E. 

l 


X 

V. 


(III # 5 # 65) 


( x 2 + »J X )* 1 + 


A 


-< X 1- oj x ) 


\ 



u. 


1 



> 

-V.J 


E. 

1 


' 1 

u. 

X 


z* 

^ v i. 


(III. 5.66) 


where X now includes the Fermi energy x . It is easy to show 

that 

ti. = \ . 

i , i 


u, = v* v, = u* i. 

i ii 


-E, 

i 


(III. 5.67) 


So only err. (III. 5. 65) need to be solved. 

18 

Goodman demonstrated that if the self-consistent 
symmetry a x is employed, the HFB quasiparticle vacuum at any 
angular velocity may be written as 


| Q > = H a + n (U„ +■ V. at a±)lo> (III. 5. 68) 

° • “e/c. 6 ess 

where 

4= {V i ' (iii - s - 69) 

3 6 = 1 V ^ 


(III. 5.70) 





go that | 5 > and | g > are eigenvectors of a with eigenvalues, 
respectively, of -i and +i. 

Ill . 5. 4 The CHFB calculations in the case of pairing— 
plus-c;, a interactions 

We have made the following approximations: 

(i) The contribution of the pairing term of the 

Hamiltonian (III. 2. 3) is neglected in constructing 
the potential r given by eq. (III. 5. 32). 

(ii) The exchange part of the q.q interaction is also 
dropped in constructing r as well as A* 

We now give explicit expressions for the matrix elements 
of r as well as a i- n the " o^-oriented" basis discussed above. 


III. 5.5 


j K. > the s 
| K> = [| j 


Matrix Elements of Hartree-Pock Hamiltonian 

Let (K > denote a state in the direct basis and 

tate in the 'con jugate* basis; that is 

3' a + "1 ] ( j _ m ) >1= ’! +> (ill .5.71) 

m >+ (-1) a |U “ ' a J 1 
a a 


+ % |(j 7 -m )>1= I - 


I K> -C MoV l - l! “ 

thorj the matrix elements of the HP potential is 


(III. 5.72) 


’ K 1' K 2 


= Z 

e 


:K. , BlV|K 


6 > 




K3/K4 


Denoting the r matrix constructed in the 


K-basis for protons 


by r +,P 


we have 



* b 


p p. • p p 


p p, ,p p 


K 3 K 4 


+ + V + +> P_,_ „ * <*+ -IV + -> P ^ 

+/P “/ P 


p n p n p n p n 

+ < + + v + + > p , +<+-.|vi+->p n 

+, n -/ n — 1 


(III. 5. 73) 


where p over + (or -) implies ‘proton and n implies ’neutron 1 


<K llo + J K 2^ = | lv e I 2 X D 8K 2 


(III. 5. 74) 


Pairing matrix element 


The pairing matrix element is 


<e i' e i 1Vp,B 2b2 > - £ D \K 1 | 2 I V 1 K 2 ^ 3 D 3 2 K 3 ^ 4 %K 4 


<K 1 K 2l V pl K 3 K 4 > 


(III. 5. 75) 


The matrix elements of the pairing interaction are given by 


< j l m l / J *2 m 2 J V P l J 3 m 3 ' J '4 m 4 >= (~G) (-1) 


Jl + J3- m l- ra 3 


5 t J 2 5j 3 j 4 6m 2'- m l 6m 4'- m 3 ' 


(III. 5.76) 


Employing eg. (III. 5.76) we obtain, after some algebraic 


manipulations 


<(3 1 ® 1 bp I ^2 & 2 > G b D 3 1 K 1 D £iKl 


(-1) f 


(l D EH £ (-1) ) 

K 3 2 K 2 e '2 K 2 


where |k 1 > is given by eg. (III. 5.71) 
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III. 6.1 CHFB results for the doubly even Molybdenum 
i sot opes in the (3s.jy 2 f 2d 2 ^ 2 ^ 2( ^5/2 / lg 7/2 ' 

^11 y"v snare ' 


In this section we present the HFB as well as 

92- 104 

the CHFB results for the nuclei Mo. 


The choice of the valence space as well as the inert 
core is the same as in our earlier CHF calculations. The 
strengths of the various corqponents of the pairing - plus - 


quadrupole-quadrupole interaction are: 

V <=W - - 0 - 016 HsV b " 2 * 4 * * * * * ' Xnc 


= -0.0352 MeV b" 


and G =0.18 MeV. 

P 

Table III. 6 summarizes the HFB results for the Mo 

isotopes, A comparison of these results with the ones given 

in the Table III. 2 brings out the role of the like-particle 

2 

pairing correlations. It turns out that the < Q q > hfb values 

2 

are only marginally different from the <0^^ values given 

earlier. The calculation errphasizes the role of the lh ^^/2 

orbit in the context of the onset of large deformations; the 

i . 100-104.. 

lhn/2 occu P at; ^ on numbers for the deformed nuclei Mo 

are 2.01, 3.86 and 4.59 respectively. The 2dy 2 occupation 
numbers, which are in fact slightly smaller than their HF 
values, further underscore the result that the occurrence 
of large deformations in this region is not compatible with 
the usual assumption of an inert Sr fe.=38, N=56) core. 

The CHFB spectra for the nuclei 10 ^Mo have been 

given in figure 4. The present calculation reproduces 
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Fig. 4. Comparison of the calculated and the experimental 
yrast spectra in the nuclei ,0A Mo and ,0S Mo. 
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TABLE III. 7 

104 • • ■ 

Details of Mo wavef unctions. The occupation probabilities 

2 

v a as well as the pairing gaps a are given for four neutrons 
pairs near the Fermi surface ?or the various yrast levels. 
Here oj denotes the cranking frequency required to obtain the 
various J's. For each J the states a are ordered with 
decreasing v^ . 


J ir 


0) 


<J X > a 



A ~(MeV) 
aa 


2 + 0.1129 2.444 


4 + 0.1789 4.467 


6 + 0.2467 6.481 


1 

0.984 

0.51 

2 

0.958 

0.54 

3 

0.607 

0.53 

4 

0.306 

0.50 

1 

0.984 

0.50 

2 

0.962 

0.53 

3 

0.622 

0.52 

4 

0.286 

0.49 

1 

0.986 

0.49 

2 

0.967 

0.52 

3 

0.652 

0.51 

4 

0.254 

0.48 
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TABLE III. 8 

106 

Details of Mo wavef unctions. For further comments, 
see Table III. 7. 


w <J > 

X 


2 + 0.1223 . 2.446 


4 + 0.1923 4.478 


6 + 0.2630 6.483 


a 

2 

v a 

t aS (MeV) 

i 

0.986 

0.45 

2 

0.915 

0.48 

3 

0.834 

0.43 

4 

0.195 

0.47 

1 

0.987 

0.45 

2 

0.928 

0.47 

3 

0.812 

0.42 

4 

0. 202 

0.46 

1 

0.988 

0.45 

2 

0.939 

0.46 

3 

0.789 

0.42 

4 

0.204 

0.45 
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104 

successfully the yrast spectrum in Mo. The computed yrast 
106 

spectrum in Mo is, however, found quite compressed compare 
to the observed one. In Tables III. 7 and III. 8 we have 
presented some of the details of the yrast levels in the nucl 
104, 106 Mq ^ a close scrutiny of the values reveals a smal 

but systematic reduction in the pairing correlations in going 
from J v = 2 + to J 1T = 6 + in 104 Mo. 

III. 6. 2 CHFB calculations for the doubly even Zirconium 
isotopes in the (2p 1 / 2 , 3s*Jy 2 , 2< ^3/2 / 2<a 5/2' 

~VV2> lg 9/2' lh ll/2> 3pac * 

We present here first the HFB results for the 

doubly even Zr isotopes. The valence space includes, in 

addition to the orbits employed in the calculations just 

discussed, also the orbital Sp^. The strength parameters 

characterizing the effective interaction are: x nn (=Xp p ) =-0.01 

MeV b” 4 , X =-0.0231 MeV b“ 4 and G = 0.18 MeV. In addition 
np p 

to the SPE's quoted in section III. 2, we have used s (2p^ 2 ^ = 
-0.8. MeV. 

In Table III. 9 we present the intrinsic quadrupole 

deformations as well as the sub— shell occupancies for the 

94-102 

intrinsic HFB states in the nuclei “ Zr. It is seen thai 
the results are not qualitatively different from the ones 
obtained from the HP picture. The calculation reproduces 
the sudden onset of deformations in going from the nucleus 
98 Zr to 10 °Zr; whereas the intrinsic HFB quadrupole moment 
in 98 Zr is only 39.3b 2 out of a maximum of 76.7 b , that 
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Zr is approximately 64 per cent of the maximum 
possible value for the given space. 

In figure 5 we have displayed the sub-shell occupation 
numbers for the various orbits that we have calculated by 
employing the HFu intrinsic states for the nuclei 94 “ 102 Zr. 

The observed spectrum (see Fig, 1) in ^Zr is characterized 
by a rather large (^1,76 MeV) value of the (E 2 + - Eq+) separa- 
tion, It is reasonable to interpret this feature in terms of. 
the 2d^y2 sub-shell closure. The computed values for the 
2 dg /2 occupation numbers are, in fact, seen to display this 
effect. 


The computed values for the 1^21/2 occupancies bring 

out the strong correlation between the increased occupation 

of this level and the sudden onset of large deformations in 

100'102 Zr ^ Further, the present results indicate that the 

onset of deformations is also marked by a significant depletior 

of the 2dy 2 orbit; the sub-shell occupation number for this 

orbit in 100 Zr is only 3,46, This invalidates the usual 

94 

assumption concerning the inertness of Sr core which has 

. . 4—6 

formed the basis of the earlier investigations regarding 
the origin of deformations in the Zirconium region. 


Federman and Pittel 5 have earlier conjectured that the 

100 . 102 „ , . . 

occurrence of deformations m Zr may be arising 

primarily from the p-g 9 / 2 )p “ ^ lg 7/2) n ] effective interaction, 
This mechanism demands a correlation between the 1 g9/2“P roton 
and lg- 7 /o -neutron occupation numbers. In view of the small , 
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corrputed values of the occupation numbers - the values 

are seen to range only between 1,36 and 2.44 in the nuclei 
9 S— • 102 

Zr - the present calculation does not support this 


idea. 


Table III. 10 summarizes our HFB results for the isotopes 

92-106 2 

Mo. The corrputed < Q values are quite consistent 

O rib id 

with the observed deformation trends (see Pig. 1(b)). The 

(E 2 + - Eq+) separations in the nuclei ^^Mo are 0.19 MeV 

and 0.17 MeV respectively; this is to be contrasted with the 

value 1.51 MeV for this separation in the nearly. spherical 
92 

nucleus Mo. The present HFB calculation yields intrinsic 

deformations for 'the nuclei -^4, 106 ^ q w hich are about 96 percen 

2 

of the respective <Q 0 > maX values. 


A consideration of the lfr^l /2 numbers again 

reveals the connection between them and the onset of large 
deformations in Mo isotopes. Further an examination of the 
lg 9/2~ proton and the lg 7/2~ neutXdn occu P ation numbers also 
fails to suggest any increased involvement of the j^lcfgy^) p 

(lg ?/2 ) n ^j effective interactions in the quasirotational 
100-106,, 

isotopes Mo. 


In figure 6 we present the CHFB yrast spectra for the 
nuclei In ^^Zr the calculated (E^+ - E 4 +) value 

is- about 0.5 MeV too large compared with the experiments. ' The 
yrast spectrum in "^^Zr is# however, in fair agreement with 


the observed one 



well as the sub-shell occupancies for the HFB intrinsic stc 


< 

H 




0 ) 

£ a 

C-i -a 


03 CM 

o LO 
4-> 'd 
O CM 
Pi 

0 . CM 
> 
CM 


U JCO 
O 


O 

G* 

m t -4 

a 

O CM 

M \ 

4-> LfT 

d 'd 

0 CM 

Si 

CM 




8. Comparison of the calculated and the experimental ytasf 
spectra in the nuclei 15e Zr end m Zr. 
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TABLE III. 11 

Details of wavef unctions. . The last column gives the 

i 2 

contribution to the angular momentum, J a = j~< J x > a a + v a * 

for various orbits a . * 


J * 

CD 

< J x > 

a 

< 

CMeV ) 

J a 




i 

0.957 

0.87 

0.11 




2 

0.797 

0.86 

10.01 

2 + 

0.1480 

2.451 

3 

0.376 

0.85 

0.74 




4 

0.297 

0.86 

1.99 




5 

0.214 

0.87 

- 0.38 

i 




1 

0.952 

0.88 

0.66 




2 

0.824 

■ 0.85 

10.25 

4 + 

0.1880 

4.765 

3 

0 . 360 

0.85 

1.02 




4 

0.267 

0.84 

• 3.04 




5 

0.227 

0.85 

- 0.42 




1 

0.954 

0.86 

0.36 




2 

0.859 

0.83 

10.33 

6 + 

0.217 

6.490 

3 

0.366 

0.85 

1.00 




4 

0.281 

0.83 

3.54 




5 

0.225 

0.85 

- 0.44 
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TABLE III. 12 


Details of 


102 


Zr wavefunctions 


J * 

6) 

<J x> 

a 

Va 

u (M2V) 

J a 




i 

0.957 

1.00 

. 1.63 




2 

0.831 

0.97 

9.94 

2 + 

0.1310 

2.43 

3 

0.578 

■ 0.98 

0.81 




4 

0.469 

0.98 

1.73 




5 

0.441 

1.00 

-0.27 




1 

0.957 

1.00 

1.07 




2 

0.858 

0.96 

10.04 

4 + 

0.1620 

4.466 

3 

0.574 

0.97 

0.83 




4 

0.485 

0.97 

2.49 




5 

0.436 

0.99 

-0.18 




1 

0.956 

0.99 

0.860 




2 

0.873 

0.95 

10.02 

6 + 

0.1965 

6.500 

3 

0. 556 

0.95 

0.02 




4 

0.512 

0.96 

3.14 




5 

0.457 

0.97 

-0.00 




1 

0.956 

0.98 

0.77 




2 

0.897 

0.93 

10.18 

8 + 

0.2320 

8.490 

3 

0.553 

0.94 

0.88 




4 

0.529 

‘ '0.94 

3.66 




5 

0.451 

0.96 

-0.04 
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In Tables III. 11 and III. 12 we present some of the 


details of the CHFB wavef unctions for the yrast levels in 

100, 102 


Zr. The non-axial orbits of lh^^ parentage can easily 


be identified here through their relatively large ('vlO.O) 


contributions towards <J„>. The occupation probabilities for 

A 


the non-axial h. .. /0 orbits show a systematic increase in 

11/2 '“ + 

going from the J 77 = 2 + to J 1T = b members of the yrast cascade 


• 4 .^ , . 100 , 102 _„ . 

xn the nuclei Zr 


III. 7 Conclusions 


The purpose of the present investigation was to 
examine, in a microscopic manner, the origin of the dramatic 


onset of large deformation in the Zirconium region. Some recen 
6 


investigations have attempted to established a link between 
the occ ur rence of large deformations in the region with the 
simultaneous occupation of neutrons and protons of the spin- 
orbit partner orbits lg ?/2 and lg g/2 . In view of the calcula- 
tions presented here, the earlier result appears to be an 


artefact of the model employed; the earlier calculations invol^; 


an inert ^Sr core with the protons confined to the orbits 

2 P 1/2 ' lg 9/2 and 2d 5/2 and the neutrons t0 the 3 s 1/2 ' 2d 3/2 

orbits. The following conclusions can be 


lg 7/ 2 and lh ll/2 


drawn from the HFB (as well as some CHFB) results obtained 


in the present work; 

(i) The involvement of the 1^1/2 or kit one the most 
inportant factors responsible for the observed dramatic 
onset of large deformations at N=60» 
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(ii) The observed large deformations in the Zirconium 

region imply a significant polarization of the Sr 

94 * 

core. In other words# the assumption of an inert Sr 
core is not consistent with the observed rotational 
features in the doubly even isotopes of Zirconium and 
Molybderfum with A>100. 

(iii) The calculations do not reveal any significant connection 
between the occurrence of deformations in this region and 
the selective occupation of the (lg^^) -proton and 

the (lg^ 2 ) *" neutron orbits. 
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CHAPTER XV 


MICROSCOPIC DESCRIPTION OF THE YRAST AND 
THE YRARE BANDS IN 10 °Pd AND' 102 Pd 

I V. 1 Introduction 

1 

Scharff-Goldhaber et al. have some. time ago 

reported the level schemes of the nuclei ^®pd (N=54 A z=46) 

102 

and Pd (N=56, Z=46) populated by heavy-ion induced reactions. 

100 91 12 

The levels in Pd were populated in the Zr { C , 3n) and 

the ^Zr (^ 2 C, 3n) reactions. The spectra in ^ 2 Pd were 

obtained through the Zr ( x C, 3n) and Zr ( C, 4n) 

reactions. 

The observed yrast as well as the yrare spectra in the 

nuclei ^ 0#1 ® 2 Pd are marked by considerable departures from 

the rotational behaviour (see Fig. 7). The yrast band in ^ 00 Pd 

(J =14) was found to deviate significantly from the "variable 
' max 

2 

moment of inertia" (VMI) predictions . The levels of the 
yrare band with J*i8 + also -displayed similar properties. The 
yrare band in 100 Pd also had J maX = 14 and an experimental 
search for the low-spin members (with J<8) of this band had 
proved inconclusive. The yrast band in ^^Pd, reaching J=14- 7 
was found to deviate only slightly from the VMI predictions. 
However, the available yrare states in Pd with J =10 , 12 
and 14 + showed significant deviations from even the VMI law. 




Fig. 7. Observed level schemes for the nuclei 100,102 Pd . 
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The earlier phenomenological attempts 2 ' 3 at understanding 
the observed levels in 100 ' 102 pa have met with limited success. 
Apart from the earlier studies in the framework of the VMI 
model, an attempt was also made by Smith and Volkov* to explain 
the observed features of the yrast and the yrare bands in the 
Palladium isotopes by invoking instability towards asymmetric 
deformations at sufficiently high (>8) angular momenta? However 
a lack of any microscopic calculation has hindered an under- 
standing of the observed spectra in terms of the underlying 
single-particle states. 

In this chapter we carry out a microscopic study of the 
yrast and the yrare bands in the nuclei l00 * 102 P< j by employing 
the variation-after-projection (VAP) formalism^ in conjunction 
with the HFB ans'atz for the axially symmetric intrinsic wave- 
functions. The choice of the VAP method was dictated by two 
considerations. Firstly, it is desirable to use a calculations 
framework which allows for the possibility of having different 
intrinsic states for each yrast level, in view of the observed 
large deviations of the yrast levels in these isotopes from 
the J(J+l)-law. Secondly it is rather easy to compute the 
intra-band E2 transition probabilities in the VAP method. In 
the present work we have again used the pairing-plus-quadrupole 
quadrupole effective interaction operating in the 3s-2d-lg-lh^ 1 
configuration space. 

The results presented in this chapter bring out the 
significant role played by the lh^.^ orbit In lending a 
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somewhat dichotomic character to the yrast spectra in the 
nuclei ^0, 102p^ it turns out that the yrast spectra can be 
interp reted in terms of two intrinsic states occurring at about 
the same energy but possessing different sub-shell occupancies 
for the lh^^ orbit. 

In section IV. 2 we discuss the method of calculation. The 
procedure for obtaining axially symmetric HFB intrinsic states 
(with K = < J > =0) is already contained in our earlier discussion 
of the CHFB method that we have presented in the preceding 
chapter. We have, therefore, confined ourselves to only a 
brief outline of the technique of exact angular momentum 
projection from HFB intrinsic states. In section IV. 3 we 
present and discuss our VAP results. Finally, section IV. 4 
contains some concluding remarks. 


I V. 2 Calculations! framework 

IV. 2.1 Angular momentum projection from HFB intrinsic states 
For axially symmetric HFB intrinsic states (xcLth 
K=0) the states )a> appearing in (III. 5. 51) can be labelled by 
their 'm' values, where m = <j >. Thus the HFB state can 

5 

be written as • 


® >= n (U 1m + V, m at at_)|o > 

*o . rm im im 

im 


(IV. 2.1) 


Here ‘i 1 labels the different deformed orbits with the same 
value of m. 

The creation operators at can be -written as 


a + = .id, . at ; at- = I (-1) 
im J im, 3 jm im j 


1+j-rtu + 

D im, j ^j-m 


(IV. 2. 2) 



Here j labels the spherical single-particle orbits lg, 


2d 5/2 ' 3s 1/2 ' 2d 3/2 ' lg 7/2 ' lh ll/2* 

The wavef unction e:xpressed in (IV. 2.1) 
the form 


9/2 ' 

can be recast into 


with 


j* >= N exp (■* E f a* at)]o> 
° 2 a e ctg a ^ 


(IV. 2.3) 


■ a B 


a a 


p 


uv. 


Here we have denoted the quantum numbers ( j m Q ) by a , and N 
is a normalization constant. 

Writing the shell-model Hamiltonian (III. 5.1) 


as 


H = 1 ‘s.M a 0 + i cm <C,B1v a I y6 > ^ a P a 8 a Y 


a 


we can write the energy of a state with angular momentum J as 
e .t - ««ol H P o„ IV/ <S ol P oJ ! o > 


TT 


= (/ h(e) d^ o ( e) sinede)/(/n(e) d^ o (0) sine de) 

o 0 

0 (IV. 2.5) 


The intensities of the various angular momenta contained in the 
intrinsic wavefunction are given by 

IT 

a J = i (2J+ 1 ) / n(0) d oo (e) Sin9de (IV. 2.6) 

0 

The overlap integrals h(g) and n(6) are given by 

h ( 6) = n(e) [Z e a 0V(l+ M )) aa + 0/4) a j y6 <a9 |v a Iy 5>. 
a 

(M/ (1+M) ) + s (1/1 +M) yp F p 5 (1/ (l+M) ) v „ f v 0 }] 

ya o p v p 

(IV. 2.7) 
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n 


( 6) = [det (1 + M(8))J 1//2 


(IV. 2.8) 


Here, 


■a 8 


(9) “ . E , 4® m. <9) \ K -h m 'e 

m * .ml n. n g g 


m a * m g a a 


M = F f ' 1 


8 

(IV. 2.9) 
(IV. 2. 10) 


Employing the angular momentum projected wavef unctions > = 


P MK ^ $ o * one can wr; ^ e 


:,P MkI Q 3 ,p MK > = \ (2J+1) vf' [ n<J uC J 2 J, 1 dJ UD (e) n <- 6) 

o M ii n -J 


V p O 

x 

x b 2 ( I e < a | Q 2 j 8 > p 3 (0)) Sin 6 de (IV. 2. 11) 

t 3 «8 t3 V 


The static quadripole moments as well as the E2 transition 
matrix elements, B(E2; jt + j+) , can easily be obtained by 
using the above relation. Here b is the oscillator length 
parameter and Q v = (16ir/5) 7 (r /b ) Y y (o). The matrix p 
appearing in eq. (IV. 2. 11) is defined as 

P = M/U+M) (IV. 2. 12) 

J 

Further, the normalizations n are given as 

# 

n j = j n(e) d^e) sinede (IV. 2. 13) 

o 

The method of performing the projection calculation is as 
follows. Employing the HFB wavefunctions we first set up 
the f matrix. Then F, M and 1/(1+M) are evaluated for 
twenty values of the gaussian quadrature points between the 
range (0, -n/2). (The axial symmetry of the intrinsic state 
permits us to change the integration limits to (0, ir/2) ) . 
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The projected energies are then computed from the expressions 
(IV. 2.5 - IV. 2. 10). 

IV. 2.2 The Variation-after-pro lection (VAP) method for 
studying the yrast spectra 

The VAP calculations have been carried out by 

employing the following procedure. We first generate self- 

consistent HFB solutions, 4(3), by carrying out HFB calculatio 

2 

with the Hamiltonian (H - 3Q 0 ). The intrinsic state for each 
yrast level is then selected by finding out the minimum of 
the projected energy Ej(B) (= < $ (B) j H | 0 (3)>/ < 4(3){P^ o | 
$(&)>) as a function of g . In other words, the intrinsic state 
for each J satisfies the following condition. 

s[>(3) l H P ool S(3)>/<4(3) |PoqI •(&) >3 = 0 (IV. 2. 14) 


IV. 3 Results and discussion 

The spherical SPE's that we have employed in the 
present calculation are (in MeV) : e(ig g ^ 2 ) - o.5, e(2d 5 ^ 2 )=5.4 
e ( 3s-, /r> ) = 6.4, e ( 2d 


l/2 y ~ k^^ 3 / 2 ) = 7,9/ e (lg- 7 / 0 ) = 8.4 and e (lh^ / 0 )= 8 . 


The strengths of the various components of the pairing-plus- 


quadrupole-quadrupole interaction were taken as: 


( = ^ p ) = -0.012 MeV b” 4 , x_ = -0.02265 MeV b ~ 4 and 


nn 


np 


Gp = 0.18 MeV. The present set of input parameters differs 


slightly from the earlier one employed in the CHFB calculation 
for the Molybdenum isotopes. We have reduced the ( ^ ri n/2~^9/2 
separation by 1.1 MeV. Further, the strengths x* s have been 
reduced by about 25 percent. The lowering of the lh.^^ 
energy is intended to mock up, at least partially, the, effects 
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due to the neglected, higher-lying single-particle 
orbits. 

In tables IV. 1 and IV. 2 we have given the details of our 
PHFB as well as VAP calculations for ±i ^Pd« 

The mi nimum- energy HFB solutions for 100 ' 102 pd are ■ 

characterized by large intrinsic quadrupole moments (about - 
2 2 

54b and 66b respectively). This is consistent with the 

r 

O 

experiments of L. Hasselgren et aJL. which yield some indica- 
tions concerning the prolate shapes in the Palladium isotopes 
with A= 10 2- 106. 

A comparison of the VAP as well as the PHFB results 

presented in the table IV. 1 reveals a significant change in 

the yrast spectrum in ' 1 ' 00 pd at J=2 - the intrinsic quadrupole 

2 2 

moments increases from 54.1 b to 67.3 b . The results for the 
sub- shell occupation numbers indicate that the increase in 
the deformation is correlated with a noticeable' increase in 
the lh^^ 2 ” neutron occupation at the cost of the neutrons in 
the lg g / 2 orl:> it. 

102 

The VAP results for the nucleus Pd (see table IV. 2) 

indicate a significant change of the intrinsic state in going 

from the yrast J=4 to J=6. Here again a change in the 

intrinsic deformation of the yrast cascade is accompanied with 

an enhanced (neutron) occupation of the Ih^^ or ^ t an< 3 a 

depletion of the lg g , 2 neutrons. The changes occurring in the 

2n 2 > values of the yrast levels in 100 ' 102 P d are roughly equi- 
^o 

valent to an increase of the effective Nilsson deformation 




*F*f ^ 
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PHFB as well as the VAP results in 1U,!:: pd. For further comments see Table IV. 1 
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In figure 8 we present a comparison of the VAP and the 

observed spectra in 100, I02p^ Here we notice that the observed 

level spacings in the lower part of the yrast cascade (with 

J=0 through 8) are reasonably well reproduced by the present 

100 

calculation for these nuclei. In Pd, however, the positions 

of the calculated yrast states with J=10 and 12 are about 1.5 MeV 

102 

too high compared with the experiments. In the nucleus Pd 
the computed position of the J=12 level shows a discrepancy 
of about 1.3 MeV. 


We have indicated in figure 8 some of the observed, high- 
spin yrare levels in 100,102p^ The figure also shows the 
positions of these yrare levels calculated by carrying out 
angular momentum projection on the intrinsic HFB state. The 
PHFB states are seen to be in poor agreement with the observed 


yrare levels. The results thus indicate that the mixing of the 
K=0 HFB intrinsic states with other high-lying (quasip article- 
excited) K=0 bands may be quite significant from the point of 


view of an accurate description of the high-spin yrare levels 

100, 102 , 


In tables IV. 3 and IV. 4 we have given the reduced transi- 
tion probabilities for electric quadrupole transitions, B(E2), 

as well as the static quadrupole moments for the yrast cascades 

. 100 , 102 -- 

m ' Pd. 

A comparison of the computed values of the electa. omagnetic 
observables with the experiments is expected to provide a 
sensitive check on the vavof unctions employed. However, a lack 




o 




Fig. 8. Calculated (TH) and experimental (EXPT) levels in 100, 102 Pd . 



100 


TABLE IV. 3 

The reduced transition probabilities for the electromagnetic E2 
transitions as well as the static quadrupole moments for the 
yrast levels in ^^Pd. The reduced transition probabilities and 
the static moments have been given in units of e .b and e.b 
respectively, where b is the oscillator parameter. Here e p (e^) 
denotes the effective charge for protons (neutrons) . 


Transition + j * 
CE(Jf)-E(J+0in KeV 


CB(E20 1 / 2 


Q(jJ) 


0 + 2 + 
(665.3) 


10.75 e p + 13.91 e n 


-6.86 e — 8.88 e 
P n 


2 + + / 1 
(750.5) 


6.51 e p + 9.53 e n 


• 8.88 e p -13.12 e n 


4 + 6 + 

(773.0) 


15.66 e p + 22.75 e n 


■ 11.23 e -16.24 e^ 
p n 


6 + 8 ^ 
(798.6) 


16.03 e + 23.28 e 

Jr 


■ 11.88 e p -17.12 e n 


8 + -*■ 10 + 
(881.4) 


16.23 e p + 23.57 e n -12.34 e p -17.70 e n 


10 + -*• 12 + 
(892.0) 


16.35 e p + 23.75 e n 


■ 12.69 e p -18.11 e n 


12 + -»- 14 + 
(945.4) 


16.39 e p + 23.35 e n 


■ 12.98 e p -18.42 e R 
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TABLE IV. 4 

The results for the reduced transition probabilities as well as 

10 2 

the static quadrupole moments for the yrast levels in Pd 


Transition 

T + 

J f 

Cb(E2)I] 1//2 



Q(Jp 


jj2 (J+)-E(jtQ in 

KeV 










0 + -v 2 + 

(556.4) 


12.19 e 

P 

+ 18.01 

e n 

- 7.79 

e -11.57 
P 

e n 

2 + -»■ 4 + 

(719. 3) 


14.57 e 
]? 

+ 21.69 

e n 

-9.94 

e -14.73 
P 

e n 

4 + -► 6 + 

(835.5) 


6.17 e p 

+ 8. 35 

e n 

-12.72 

e -18.16 

ir 

e n 

6 + 8 + 

(901.4) 


19.09 e 

ir 

+ 27.76 

e n 

-14. 19 

e p -20.29 

e n 

8 + + 10 + 

(714.6) 


19.19 e 

+ 28.09 

e n 

-14.65 

e p -20. 90 

e n 

10 + * 12 + 

(704.8) 


19.17 e 

ir 

+ 28 . 28 

e n 

-14.96 

e -21. 36 
P 

e n 

12 + 14 + 

(892,4) 


19.07 e D 

+ 28.38 

e n 

-15.19 

e —21.61 
P 

e n 
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of availability of the yrast life-times in "'■^ 2 pd has 

prevented us from testing the VAP wavef unctions obtained in 
the. present calculation. 

We have noted earlier the significant increase in the 
deformations of the intrinsic states at the J=4 and J =6 yrast 
levels in •'•^pd respectively. From the B(E2) results 

presented in the tables IV. 3 and IV. 4 we find that these 
changes in the intrinsic deformations manifest themselves here 
in terms of sizable reductions in the B(E2) (2 + -*'4 + ) and 
B(E2) (4 + -* 6 + ) values in • ! - 00 '^ ( - ) ^pd respectively. 

100 102 

A measurement of the yrast half-lives in * Pd would 
not only be interesting in itself, it would also help establish 
the band parentage of the various yrast levels in these 
nuclei. 

IV. 4 Conclusions 

Summing up, the present VAP calculations for the 
low-lying levels in - L °°'. 102 Pd are fairly successful in explain- 
ing, in a microscopic manner, the significant departures from 
the rotational behaviour that characterize the observed yrast 
spectra in these nuclei. 

The results presented here reveal that the increase in the 
intrinsic deformations of the yrast levels (with J>_2) in 
100,102 pd is due tQ the enhanced occupation of the 
orbit. In other words, the observed "centrifugal stretching" 
in 100,102 pd is due tQ a selective participation of the lh 13 / 2 

orbit for some of the yrast levels. 



103 


REFERENCES 

1* G. Scharf f-Goldhaber, M. McKeown, A. H. Lunpkin and 
W.F. Piel Jr., Phys. Letters 44B (1973) 416. 

2. M.A. J. Mari sco tti, G, Scharf f-Goldhaber and B. Buck, 

Phys. Rev. 178 (1969) 1864. 

3. B.C. Smith and A. B. Volkov, Phys, Letters 47B (1973) 193. 

4. L. Satpathy and S.C.K. Nair, Phys. Letters 2613 (1968) 257, 
R. Dreizler, P. Federman, B. Giraud and E. Osnes, Nucl. 
Phys. All 3 (1968) 145; 

J. Vojtik, Nucl. Phys. A212 (1973) 138; 

E. Caurier and B. Grammaticos, Nucl. Phys. A279 (1977) 333. 

5. N. Onishi and S. Yoshida, Nucl. Phys. 80 (1966) 367. 

6. L. Hasselgren, C. Fahlander, J.E. Thun, B. Orre and 
N.G. Jonsson, UUIP-957 (1977); 

L. Hasselgren, in Interacting Bosons in Nuclear Physics ,- 
edited by F, lachello (Plenum Press* New York and London, 
1979), p 67. 



104 


CHAPTER V 
CONCLUSIONS 


Looking through the chart of the nuclides one comes across 
an interesting region between' A=90 and A=110. In this region 

one encounters at one end classic shell— model nuclei such as 

90,, , 92 102 

Zr and Mo and, at the other end, classic rotors like Zr, 

104 , 106 

Mo and Mo. The present study was motivated by a desire 
to gain some insight into this phase-transition from spherical 
to deformed nuclei. 

Our HP-oriented study of the effective interactions 
operative in the 2s-ld, 2p-lf and 2p-lf-lg shells provided a 
qualitative understanding of the occurrence of deformations in 
the Zirconium region. An important spin-off of this study was 
the result that the relative importance of the pairing correla- 
tions in the like-particle and the neutron-proton channels 
depends very much on the spectrum of the spherical single- 
particle energies. An appreciation of this role of the single- 
particle energies turns out to be instrumental in understanding 
the absence of rotational features in the N>Z nuclei in the 
2s- Id and 2p-lf shells. 

In chapters III and IV we have focussed our attention on 
the quantitative description of the qua si rotational yrast 
spectra in some doubly-even Zirconium, Molybdenum and Palladium 
isotopes, enploying the CHFB and the VAP techniques in 
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conjunction with the pairing-plus- quadrupole- quadrupole effect- 
ive interaction. The calculations indicate that the origin of 
the sudden onset of large deformations in the highly neutron- 
rich isotopes of the zirconium and Molybdenum can be traced 
back to the involvement of the 1^21/2 or ^ t the underlying 
configuration space. The observed centrifugal stretching in 
the nuclei 100, 102p£ j_ s a j_ so found to be connected with an 
enhancement of the ^ 2 . 1/2 occvl P a -' t: ^ on (for neutrons) in the 
J>2 states of the yrast bands. 

The present study was concerned primarily with the rota- 
tional features in the available yrast spectra with J<j3 in the 
Zirconium region; we did not examine the high-spin spectra in 
this region. The calculations presented here have revealed 
that the time-reversed pairs of 1 ^ 2 . 1/2 neutrons are just below 
the relevant Fermi surfaces in the Zirconium .and Molybdenum 
isotopes. Therefore, at higher excitation energies, the orbit 
lhn /2 is likely to play a role somewhat akin to that played 
by li^3/2 or kit i- n t ' iie ^are-earth region where it leads to 
backbending effect through a realignment of the 11^3/2 P a ^ r * 

It will thus be very interesting to study the high— spin yrast 
spectra in the Zirconium region, both from theoretical as well 
as experimental point of view. 
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